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1.0  INTRODUCTION 
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1.1  Development  of  the  PIS  Ratio  Discriminant 

The  regional  PIS  amplitude  ratio  has  long  been  investigated  as  a  possible  discriminant 
between  explosions  and  earthquakes,  because  intuitively,  earthquakes  being  shear,  dislocation  type 
sources  should  produce  more  shear-wave  energy  than  explosions  (Pomeroy  et  al,  1981).  Thus, 
PIS  ratios  should  be  lower  for  earthquakes  than  for  explosions.  Early  studies  in  the  1970’s  and 
early  1980's  demonstrated  the  potential  of  the  discriminant  (e.g.,  Blandford  et  al,  1981)  for 
separating  explosions  and  earthquakes.  However,  studies  of  the  PILg  ratio  discriminant  at  low 
frequencies  by  Nuttli  (1981),  Bennett  and  Murphy  (1986),  and  Taylor  et  al  (1989)  found 
significant  overlap  in  the  earthquake  and  explosion  populations  for  events  in  Eurasia  and  in 
Western  United  States. 

Bennett  et  al  (1989)  demonstrated  improved  separation  of  nuclear  explosions  and 
earthquakes  on  the  basis  of  high  frequency  SIP  ratios.  Dysart  and  Pulli  (1990)  showed  the  same 
trend  for  mine  blasts  and  explosions  in  Scandinavia  recorded  at  the  NORESS  array.  However,  in 
both  studies,  the  earthquake  and  explosion  groups  were  not  in  the  same  geographic  region  and 
propagation-path  differences  may  have  biased  the  discrimination  results.  Baumgaidt  and  Young 
(1990)  studied  mine  blasts  and  earthquakes  which  were  in  the  same  region  in  Scandinavia  and 
showed  that  mine  blasts  have  significantly  higher  PnlLg  ratios  at  high  frequency  than  earthquakes. 
Interestingly,  they  showed  that  the  discrimination  was  due  to  explosions  having  higher  Pn  energy 
than  to  earthquakes,  whereas  their  Lg  levels  were  about  the  same.  Kim  et  al  (1993)  also 
demonstrated  the  success  of  this  discriminant  at  high  frequencies  for  quarry  blasts  and  earthquakes 
in  Eastern  U.S. 

1.2  The  Discriminant  Transportability  Problem 

Discrimination  research  has  always  been  burdened  by  the  need  to  understand  the  effects  of 
propagation  path  structure  on  high-frequency  seismic  waveforms.  To  remove  the  effects  of 
propagation  path  bias,  earthquakes  and  explosions  in  the  same  geographic  region  have  been  sought 
(e.g.,  Murphy  and  Bennett,  1982;  Bennett  and  Murphy,  1988;  Baumgardt  and  Young,  1990).  In 
other  studies,  the  various  kinds  of  sources  (nuclear  explosions,  earthquakes,  mine  blasts)  have 
been  in  entirely  different  regions  (Taylor  et  al,  1988,  1989;  Pulli  and  Dysart,  1990)  which  leaves 
open  the  question  of  how  much  differences  in  propagation-paths  biased  the  discrimination  results. 
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Lynnes  and  Baumstark  (1991),  for  example,  demonstrated  propagation  path  effects  on  PIS 
amplitude  ratio  discriminants  and  Lg  spectral  ratio  discriminants  for  numerous  events  in  North 
America  that  caused  poor  discrimination  results.  One  may  also  encounter  positive  but  false 
discrimination  results  because  of  effects  of  propagation  path  if  explosions  and  earthquakes  occur  in 
different  regions  which  have  different  propagation  path  effects. 

A  good  example  of  how  propagation  path  effects  may  bias  discriminants  was  the  recent 
experiment  to  identify  a  seismic  event  near  Novaya  Zemlya,  which  occurred  on  December  31, 
1992  (Ryall,  1993).  Initially,  all  of  the  investigators  who  participated,  including  ENSCO 
(Baumgardt,  1993b),  identified  the  event  as  mine  blast,  based  mainly  on  analysis  of  high- 
frequency  Pn/Sn.  Comparison  of  the  measured  ratios  with  reference  events  placed  the  31 
December  1992  event  in  the  Kola  mine  blast  population.  Fisk  et  al  (1993)  determined  that  this 
conclusion  was  statistically  valid,  based  on  analysis  of  the  same  features.  However,  the  31 
December  1992  event  was  in  a  different  region  than  the  Kola  blasts  and  was  much  farther  away 
from  the  ARCESS  array.  Baumgardt  (1993b)  concluded  that  this  event  was  an  earthquake,  based 
on  these  noted  propagation  path  differences  and  a  relative  scaling  argument  using  an  earlier 
earthquake  recorded  at  NORESS. 

This  example  may  be  typical  of  a  future  CTBT  or  NPT  monitoring  problem,  where  an 
event  of  unknown  identity  occurs  in  a  region  where  there  has  been  little  or  no  seismicity  or  a 
limited  seismological  record  available.  To  try  to  characterize  such  events,  which  we  call  "special 
events,"  feature  comparisons  must  be  made  between  the  unknown  event  and  known  events  in  other 
parts  of  the  world.  This  "discriminant  transportability"  problem  requires  the  monitoring  system  to 
account  for  such  propagation  path  biasing  effects,  to  correct  for  them,  or  to  use  discriminants 
which  are  insensitive  to  such  biasing  effects. 

There  are  two  problems  associated  with  transporting  the  P/S  ratio  discriminants  (Pn/Sn, 
Pn/Lg)  from  one  region  to  another. 

(1)  These  amplitude  ratios  depend  on  distance  of  propagation,  even  in  the  same 
tectonic  region,  because  the  Pn  and  Lg  amplitudes  attenuate  at  different  rates 
with  distance.  The  same  is  also  true  of  the  Pn/Sn  ratio,  used  in  regions 
where  Lg  waves  are  blocked,  such  as  for  events  at  Novaya  Zemlya  recorded 
at  ARCESS  (Baumgardt,  1993).  Thus,  classifying  seismic  events  with 
reference  or  training  events  at  different  distances  will  require  some  kind  of 
path  correction  for  ^stance  which  accounts  for  the  differential  attenuation  of 
Pn  and  Sn  or  Lg. 

(2)  Pn/Lg  ratios  can  vary  significantly  due  to  partial  or  complete  blockages  of 
Lg,  which  Baumgardt  (1990)  showed  could  be  caused  by  sudden  variations 
in  shallow  crustal  structure,  such  as  when  Lg  waves  propagate  across 
sedimentary  basins.  Lg  "blockages"  can  cause  sharp  increases  in  the  Pn/Lg 


2 


ratio  and  could  make  earthquakes  have  "explosion-like"  high  ratios  which 
could  cause  them  to  be  misclassified  if  the  blockages  are  not  identified. 

This  report  presents  the  results  of  a  study  which  addresses  these  two  problems  associated 
with  transporting  regional  discriminants.  First,  to  transport  PIS  amplitude  ratios  from  one  distance 
to  another  in  the  same  tectonic  region,  we  present  a  technique  of  using  station  pairs  at  different 
separation  distances  to  determine  amplitude  ratios  as  a  function  of  distance.  The  second  problem 
of  correcting  for  possible  blockage  effects  in  the  upper  crust  when  transporting  discriminants 
across  tectonic  boundaries  requires  knowledge  of  the  crustal  structure.  To  correct  for  this  kind  of 
blockage,  we  have  investigated  correlations  between  PnlLg  amplitude  ratios  with  crustal 
parameters  derived  from  GIS  databases  (Fielding  et  al,  1992)  using  techniques  similar  to  those  of 
Zhang  and  Lay  (1993a,b),  and  Zhang  et  al  (1994). 

Section  2.0  of  this  report  describes  the  method  for  measuring  regional  PIS  ratios,  used  in 
the  data  analysis  for  this  study.  Section  3.0  discusses  two  methods  for  correcting  regional  PIS 
ratios  for  distance,  one  based  on  a  regional -phase  detection  study  of  Sereno  (1991)  and  a  new 
method  we  have  developed  applicable  to  regions  where  station  coverage  is  sparse.  Section  4.0 
addresses  the  more  difficult  problem  of  transporting  the  PnlLg  ratio  discriminant  to  a  new  region 
for  which  there  is  limited  historical  seismic  data.  Following  on  the  work  of  Baumgardt  (1990)  and 
Zhang  et  al  (1994),  we  investigate  a  method  of  using  crustal  cross  sections  between  a  source  and 
station  to  calibrate  PnlLg  ratios.  Section  5.0  summarizes  the  conclusions  and  makes 
recommendations  for  an  approach  for  transporting  the  regional  PIS  ratio  discriminant  from  one 
region  to  another. 


2.0  PIS  MEASUREMENT  APPROACH 


Following  on  the  research  of  Baumgardt  and  Young  (1990),  ENSCO  implemented  the 
regional  PIS  ratio  discriminant  in  the  Intelligent  Seismic  Event  Identification  System  (ISEIS) 
(Baumgardt  et  al.  1991).  The  system  was  specifically  designed  to  exploit  multichannel  high- 
frequency  data  recorded  at  the  regional  arrays  NORESS,  ARCESS,  FINESA,  and  GERESS.  For 
this  study,  we  rely  on  measurements  of  PIS  ratios  made  in  the  ISEIS  system  on  high-frequency 
data  from  both  arrays  and  single  stations. 


ISEIS  measures  several  varieties  of  PIS  ratios  on  incoherent  beams  computed  from 
individual  traces  recorded  at  the  regional  arrays.  An  example  of  a  bandpass  filtered  trace  recorded 
at  the  array  ARCESS  and  its  corresponding  log-RMS  incoherent  beam  are  shown  in  Figure  1. 
Given  a  set  of  seismogram  time-series  amplitudes  at  time  i  on  array  channel  Aj  (ti,  4^,  which  have 
been  prefiltered  in  the  band  Af,  the  log-RMS  incoherent  beam  value  at  time  r  is  computed  as 
follows: 


AX<.r.&f)  =  ^'Zlog 

^  y=i 


I 


(1) 


where  K  is  the  number  of  time  points  in  time  window  from  tl  to  t2,  N  is  the  number  of  channels  in 
the  array,  and  T=  These  measures  are  made  for  adjacent  time  windows  shifted 

down  the  trace,  and  when  plotted  as  a  function  of  the  window  time,  T,  the  result  is  an  envelope  of 
the  trace,  as  shown  in  Figure  1.  We  have  experimented  with  a  number  of  different  kinds  of 
incoherent  beams  and  have  found  that  the  log  RMS  display  gives  the  most  informative  view  of  the 
regional  waveform  shape.  By  averaging  log  RMS  amplitude  estimates  over  the  array,  local 
receiver  function  effects  at  the  individual  array  sites  are  averaged  out  and  a  smoothed  estimate  of 
the  seismogram  shape  results.  In  the  case  of  single-channel,  or  three-component  data,  the  same 
time-average  calculations  can  be  done  for  the  single  channels. 

ISEIS  measures  these  incoherent  beams  for  seismograms  prefiltered  in  several  frequency 
bands  using  6th  order  Butterworth  recursive  bandpass  filters.  For  high  frequency  data  at  the 
NORESS-type  arrays,  which  are  sampled  at  40  cps,  the  filter  bands  are  0.5-2.0  Hz  (for  Rg 
detection),  2-4  Hz,  2.5-4.5  Hz,  3-5  Hz,  4-6  Hz,  5-7  Hz,  6-8  Hz,  8-10  Hz,  and  8-16  Hz.  These 
pre-filters  can  be  modified  for  other  sampling  rates.  For  example,  to  process  three-component 
Chinese  Digital  Seismic  Network  (CDSN)  on  the  midband  channels  (bz,  bn,  be),  filters  designed 
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bandpass  tillered  in  me  8- lb  Hz  band.  1  ime  picks  made  by  the  analyst  are  shown,  lime  traces  are 
bandpass  filtered  in  9  bands,  (b)  Envelope  traces  are  computed  from  the  filtered  seismogram  by 
averaging  the  rms  amplitudes  in  1  second  windows  shifted  down  the  trace.  For  arrays,  IBEAMs  are 
computed  from  averages  of  the  envelopes  across  the  array.  For  single  stations,  IBEAMS  are  the 
envelopes  of  the  single  traces. 


for  20  cps  sampling  rates  are  provided  but  without  the  8-16  Hz  filter  since  this  band  exceeds  the 
Nyquist  frequency  of  10  Hz. 

Figure  2  shows  an  example  of  the  entire  suite  of  bandlimited  incoherent  beams  computed 
for  an  ARCESS  recording  of  a  Kola  Peninsula  mine  blast.  Each  of  the  incoherent  beams  have 
been  shifted  for  viewing  purposes.  These  displays  have  been  very  useful  for  analyzing  the 
frequency  trends  of  regional  seismogram  shapes  and  for  characterizing  various  regional  phases  as  a 
function  of  frequency. 

After  computing  the  incoherent  beams,  the  next  step  is  to  estimate  the  regional  phase 
amplitudes  off  the  incoherent  beams.  Regional  phase  picks,  as  shown  in  Figure  1,  are  noted  on 
the  incoherent  beams.  These  phase  picks  generally  come  from  measurements  made  in  the 
Intelligent  Monitoring  System  (Bache  et  al,  1990)  which  are  later  reviewed  by  the  analyst  We  also 
review  the  phase  picks  in  the  ISEIS  system  and  make  corrections,  if  necessary.  Amplitude 
measures  are  then  made  in  time  windows  on  the  incoherent  beams  beginning  at  these  phase  picks. 
The  time  windows  are  either  set  automatically,  using  assumed  group  velocities,  or  interactively  by 
the  analyst.  An  example  of  the  measurement  of  the  amplitude  of  a  Pn  phase  in  the  8  to  10  Hz  band 
is  shown  in  Figure  3.  Two  kinds  of  amplitudes  are  measured:  Amax,  which  is  the  maximum 
amplitude  in  the  incoherent  beam  time  window,  and  Aavg,  which  is  the  average  time  in  the 
window.  These  estimates  are  made  for  all  the  associated  phases  on  all  bandlimited  incoherent 
beams  and  are  stored  in  the  Oracle  database. 


Finally,  amplitude  ratio  measures  for  a  regional  P  and  S  phase  are  made  on  all  unique 
combinations  in  each  frequency  band: 


y?-*(/)  = 


<(/) 


(2) 


for  the  average  amplitude  estimates  and 


/?" 


....  a:„(/) 
KM) 


(3) 


for  the  maximum  amplitude  estimates.  Note  that  ratios  are  made  in  the  same  frequency  bands.  All 
these  ratios  are  computed  for  all  unique  combinations  of  phases  in  all  frequency  bands  available 
and  the  results  stored  in  the  Oracle  database. 
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times  on  the  incoherent  beams. 


Figure  3:  (a)  Amplitudes  are  measured  from  incoherent  beams  on  “phase  selections”,  or  time  windows 
of  specified  length  (AO,  starting  at  the  time  of  the  picked  phase  onset,  (b)  Both  maximum  and  average 
estimates  of  IBEAM  amplitude  are  made.  Maximum  amplitudes  are  usually  best  since  they  are  less 
dependent  on  assumed  window  length,  N. 


Measurements  of  this  kind  have  been  made  for  many  examples  of  mine  explosions,  nuclear 
explosions,  earthquakes,  and  quarry  blasts.  A  variety  of  scatter  plot  displays  are  provided  for 
displaying  the  results  in  comparison  with  the  CURRENT  EVENT,  which  is  a  user  selected  special 
event  which  needs  to  be  characterized.  Amplitude  ratios  can  be  plotted  as  a  function  of  geographic 
region,  distance,  frequency,  magnitude,  depth  of  focus  (if  known),  and  signal-to-noise  ratio.  For 
the  latter  display,  the  same  amplitude  measures  are  made  on  background  noise  and  signal-to-noise 
ratios  are  computed  for  all  phases  and  stored  to  the  database.  Examples  of  these  kinds  of  plots  will 
be  shown  later  in  this  report. 
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3.0  DISTANCE  DEPENDENCE  OF  THE  PIS  RATIO 

DISCRIMINANT 


Previous  studies  of  amplitude  changes  of  regional  phases  as  a  function  of  distance  (e.g., 
Chun  et  al,  1987;  Chun,  et  al,  1989;  Sereno,  1991)  have  clearly  shown  that  the  amplitudes  of 
regional  P  and  S  phases  will  not  decay  at  the  same  rate.  Most  of  these  studies  have  assumed  that 
the  amplitudes  of  individual  regional  phases  decay  with  distance  due  to  anelastic  attenuation  and 
geometric  spreading,  both  of  which  will  be  different  for  P  and  S.  Thus,  we  might  expect  that 
curves  of  PIS  amplitude  ratios  versus  distance,  for  events  in  regions  where  anelastic  attenuation 
and  geometric  spreading  effects  are  laterally  homogeneous,  to  be  relatively  smooth  and  systematic. 
However,  Kennett  (1993)  studied  the  distance  dependence  of  Pn,  Sn,  and  Lg  amplitudes  from  a 
single  explosion  source,  using  data  from  the  Fennolora  long-range  refraction  study  in  Sweden,  and 
found  significant  scatter  in  Pn/Sn  and  Pn/Lg  amplitude  ratios  as  a  function  of  distance.  He  argued 
that  the  complexities  of  Pn  and  Sn  interactions  in  the  crust  and  upper  mantle  cause  complex 
variations  in  amplitudes  of  these  phases  as  a  function  of  distance.  Thus,  PntSn  and  PnlLg 
amplitude  ratios  may  not  be  predictable  from  simple  parametric  expressions  for  anelastic 
attenuation  and  geometric  spreading. 

In  this  section,  we  address  the  question  of  correcting  Pn/Sn  and  Pn/Lg  amplitude  ratios  for 
distance  by  testing  the  hypothesis  that  simple,  parametric  correction  curves  can  be  utilized. 
Specifically,  we  compare  two  different  methods  for  correcting  amplitude  ratios  for  distance- 
dependent  attenuation.  The  first  one,  based  on  the  multiple-event  amplitude  distance  corrections 
determined  by  Sereno  (1991),  was  recently  used  by  Fisk  (1994)  in  a  reinterpretation  of  the  Pn/Sn 
amplitude  ratio  discriminant  applied  to  the  December  31, 1992  Novaya  Zemlya  event  The  second 
method  is  a  new  approach  developed  specifically  for  regions  where  there  are  only  a  few  stations 
separated  by  large  distances.  In  this  approach,  we  derive  a  new  set  of  Pn/Sn  and  Pn/Lg  amplitude- 
ratio  distance-correction  curves  for  Scandinavia  by  directly  fitting  the  amplitude-ratio 
measurements  themselves,  rather  than  the  individual  amplitudes,  to  simple  distance-dependent 
relations.  We  then  show  how  well  these  curves  can  correct  amplitude  ratios  for  distance.  Note 
that  both  these  approaches  assume  that  the  measurements  are  made  for  events  at  different  distances 
in  the  same  geographic-tectonic  region. 

3.1  The  Sereno  Multiple  Event  Distance  Correction  Method 

Recently,  Fisk  (1994)  utilized  the  amplitude-distance  parameters  of  Sereno  (1991)  to 
correct  Pn/Sn  amplitude  ratios  for  distance.  The  method  apparently  resolved  the  ambiguity  of  the 
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PnISn  ratios  for  the  December  31,  1992  Novaya  Zemlya  event  in  that,  after  correction,  the  ratios 
were  more  consistent  with  earthquakes.  However,  Sereno's  original  study  was  concerned  with 
detectability  of  individual  regional  phases,  not  with  distance  corrections  of  amplitude  ratios.  He 
derived  relations  for  the  distance  decay  in  different  frequency  bands  of  the  regional  phases  Pn,  Pg, 
Sn,  and  Lg  using  single-array  measurements  from  multiple  sources  of  different  magnitudes.  The 
method  involves  first  correcting  all  the  sources  which  have  different  magnitudes  to  the  same 
magnitude  by  means  of  a  frequency-dependent  source  scaling.  The  resulting  amplitudes  are  then 
fit  by  regression  to  a  set  of  analytical  expressions.  These  expressions  can  then  be  used  to  obtain 
expressions  for  amplitude  ratios  between  two  regional  seismic  phases.  Phi  and  Ph2,  as  a  function 
of  the  bandpass-filter  center  frequency,/,  in  the  form 


Ph\ 

Phi 


-(aoi-aoj)/ 


(4) 


where  n  =  af+b  accounts  for  geometric  spreading  and  the  exponential  term  accounts  for 
attenuation  relative  to  a  reference  distance,  Aq.  The  coefficients  a,  b,  and  ao  are  based  on  the 
individual-phase  distance  corrections  inferred  by  Sereno  (1991),  which  are  shown  in  Table  1. 


TABLE  1: 

Amplitude  Ratio  Distance  Correction  Parameters 
Based  on  Individual  Phase  Amplitude  Relations  of  Sereno  (1991) 


Phase 

ao 

a 

b 

Pn 

-0.02 

0.072 

1.40 

Pg 

0.11 

0.122 

2.33 

Sn 

-0.03 

0.125 

1.62 

0.19 

0.63 

0.29 
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This  formxilation  is  useful  because  it  allows  amplitude-ratio  coirections  to  be  made 
between  any  phase  pair  from  one  distance  to  another.  Fisk  (1994)  used  the  relationship  (4) 
to  obtain  distance  correction  for  the  PnISn  ratios  measured  by  ISEIS  to  a  standard  distance, 
Ao,  which  is 


^  ^  >^-0.053/-0.22 


.Sn/ 


(5) 


The  corresponding  relation  for  Pn/Lg  ratios  is 


{El 


(6) 


One  of  the  problems  of  the  Sereno  (1991)  method  of  computing  the  distance- 
dependence  of  individual  phase  amplitudes  is  the  use  of  multiple  sources  recorded  at  one 
station.  Given  the  high  variability  of  P  vs.  S  excitation  by  various  events,  the  estimation 
procedure  for  PIS  -ratio  distance  dependence  must  be  insensitive  to  variations  in  relative 
excitation.  Because  the  sources  have  different  magnitudes,  a  local-magnitude  (Ml)  source 
scaling  relationship  must  be  assumed  to  reduce  all  events  to  a  common  magnitude. 


Another  possible  problem  is  that  the  two  relations  (5)  and  (6)  above  imply  a 
completely  different  distance  dependence  of  the  amplitude  ratios,  depending  on  frequency. 
These  formulas  show  that,  for  a  given  frequency,  the  amplitude  ratios  scale  as  a  constant 
power  of  distance  relative  to  a  reference  distance.  Assuming  the  reference  distance  and 
ratio  to  be  1,  the  distance  dependences  for  3  different  frequencies,  3,  5,  and  9  Hz,  are 


f  =  3Hz  f  =  5Hz  f  =  9Hz 


Pn' 

_  ^0.379 

Pn 

_  ^0.485 

Pn~\  _  ^0.697 

.Sn_ 

.Sn. 

J 

For  frequencies  less  than  14  Hz,  the  amplitude-ratio  (R)  vs.  distance  trend  is 
predicted  by  these  relations  to  be  concave  down,  relative  to  a  slope  1  line,  as  shown  below: 


12 


As  the  frequency  approaches  14  Hz,  the  trend  becomes  a  line  of  slope  1.  Above  14  Hz,  the 
trend  is  concave  up,  as  shown  below: 


For  PnlLg  ratios,  the  trends  at  these  three  frequencies  are 


f  =  ZHz 


f  =  5Hz 


El 

Lg 


^^0.564 


El 

Lg 


f  =  9Hz 


El 

Lg 


_  ^3.912 » 


Thus,  below  3  Hz,  the  predicted  trend  for  Pn/Lg  ratios  is  concave  down,  whereas  above  3  Hz,  the 
trends  will  be  concave  up.  Thus,  for  the  main  frequency  band  of  interest  in  our  discrimination 
work,  which  is  between  6  and  10  Hz,  the  Sereno  (1990)  trends  with  distance  will  be  concave 
down  for  Pn/Sn  ratios  and  concave  up  for  Pn/Lg  ratios.  Also,  at  a  frequency  of  about  9  Hz,  the 
PnlLg  ratio  is  predicted  to  have  a  cubic  trend.  The  differences  in  the  trends  for  the  different 
frequencies  seem  extreme. 


3.2  The  Exponential  Station  Pair  Distance  Correction  Method 

A  better  method  of  obtaining  distance  corrections  for  regional  phase  amplitudes,  which 
would  not  require  assumptions  of  source  magnitude  scaling,  would  be  to  compare  the  P  and  S 
amplitudes  (or  their  ratios)  of  the  same  events  recorded  at  stations  at  multiple  distances.  However, 
because  regional  events  are  small,  most  of  them  are  usually  only  recorded  at  two  stations.  Also,  in 
many  regions,  station  densities  are  small  and  the  stations  may  be  widely  separated.  For  example, 
there  are  only  three  regional  arrays  available  for  all  of  Scandinavia,  NORESS,  ARCESS,  and 
FINESA,  which  are  each  separated  by  over  400  km.  This  may  be  typical  of  a  CTBT  monitoring 
scenario  where  only  a  small  number  of  alpha  stations  may  be  available  for  a  large  geographic 
region.  Without  using  a  smaller  density  network  of  beta  and  gamma  stations,  which  may  be  of 
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lower  quality  and  for  which  data  may  not  always  be  available,  a  method  is  needed  to  use  the  small 
number  of  stations  to  infer  distance  corrections  for  amplitude  ratios. 


We  have  developed  a  method  to  directly  infer  distance  corrections  for  frequency-dependent 
amplitude  ratios  by  measuring  the  ratios  for  events  recorded  at  pairs  of  stations  which  are  separated 
by  a  range  of  distances.  This  method  has  some  similarity  to  the  station-pair  method  originally 
developed  by  Chun  et  al,  (1987)  for  determining  the  anelastic  attenuation  of  Lg.  However, 
because  we  desire  corrections  for  amplitude  ratios,  rather  than  absolute  amplitudes,  the  reversed 
two-station  method,  developed  by  Chun  et  al  (1987)  to  avoid  biases  caused  by  differences  in 
station  calibration,  need  not  be  applied.  Thus,  we  only  require  one  event  per  station  pair,  rather 
than  two  events  needed  in  the  reversed  two-station  method.  The  amplitude-ratio  distance 
corrections  are  computed  by  a  least  squares  fitting  scheme  to  two  or  more  stations.  Thus,  for  a 
single  event  recorded  at  two  stations,  we  represent  the  amplitudes  of  the  P  (PI,  P2)  and  5(57, 52) 
at  the  two  stations  as  a  simple  attenuation-plus-geometric-spreading  model: 


P,-A  Ar"e 
1  o  1 


TtfA^ 

-n. 


TlfA. 


P^=  A  A-’^e 
2  O  ^ 


(7) 


for  the  P  waves  and 


nfA 


1 


Q  u 

S,=A  Ar^e  ^  ^ 
1  o  ^ 

TtfA^ 


5.,  =A  A 


-n. 


o 


(8) 


for  the  S  waves  at  the  two  stations.  In  these  expressions,  Ao  is  the  initial  amplitudes,  Al  and  A2 
the  distances  of  the  event  from  the  two  stations,  Qp  and  Qs  are  the  attenuation  quality  factors  for 
the  P  and  5  waves,  respectively.  Up  and  Us  are  the  group  velocities  of  the  P  and  5  waves, 
respectively,  rie  is  the  assumed  geometric  spreading  exponent  (assumed  to  be  the  same  for  P  and 
5),  and /is  the  frequency  of  the  waves.  For  the  purposes  of  our  analyses,  we  assume  that  the 
frequency /is  the  center  frequency  of  the  9  filter  bands  used  to  compute  the  incoherent  beams  on 
which  the  amplitudes  of  the  P  and  5  waves  are  measured. 


14 


Figure  4(a)  shows  a  map  of  the  locations  of  the  events  which  we  used  for  this  analysis,  that 
were  recorded  at  two  more  of  four  regional  arrays,  NORESS  (NRAO),  ARCESS  (ARAO), 
FINESSA  (FIAO),  and  APATITY  (APAO).  The  X’s  mark  the  locations  of  the  seismic  events 
which  were  analyzed.  For  each  of  these  events,  two  regional  PIS  ratios  were  measured,  PnISn  and 
PnlLg  in  the  9  frequency  bands.  The  distances  of  the  events  from  the  different  stations  ranged 
from  200  to  2000  km. 


For  station  pairs  which  recorded  single  events,  we  next  compute  the  differences  in  the 
logarithms  of  the  ratios  at  two  or  more  stations  as  follows: 


/?!  =  log  Pi  -  log5i  =  iWfAi 

=  logPj  ~  log-S^a  =  ^2 


1 


1 


QsUs 
1  1 


QsU. 


logs 

loge 


(9) 


where  R1  and  R2  are  the  log  ratios  measured  at  two  stations.  The  difference  in  ratios  is  then: 


dR  =  R2-Ri  =  ^loge 


1 

QsU. 


[^2  ^l]' 


(10) 


Because  we  have  assumed  the  spreading  to  be  the  same  for  P  and  S  waves  and  the  same  initial 
amplitude,  we  have  a  linear  relation  at  a  given  center  frequency  /between  the  difference  in  ratios 

and  the  difference  in  distances  with  zero  intercept.  The  gradient  of  the  differential  ratios, 
dR 

Gradient  =  V  =  — ,  is  then  integrated  over  specified  distance  intervals,  dA,  as  follows: 
dA 


A(A  J  = 

0 


(11) 


In  our  analyses,  we  compute  this  integral  numerically  using  distance  intervals,  dA,  of  10  km.  The 
amplitude  ratio  at  a  given  distance  is  then. 


(12) 


where  the  Slope  is  written  as 
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Regressions  To  Gradient  Integrations 
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Figure  4;(a)  Locations  of  events  and  arrays  used  to  derive  frequency  dependent  Pn/Sn  and  Pn/Lg  correction,  (b) 
Examples  of  the  distance  dependent  gradient  integration  regressions  in  the  8-10  Hz  frequency  band  for  Pn/Sn  and 
Pn/Lg  measured  for  the  events  shown  in  (a). 


Slope  =  Tif  \oge 


1 


Q  u 

PI 


(13) 


Figure  4  (b)  shows  examples  of  this  integrated  log  amplitude  ratio  as  a  function  of  distance 
for  Pn/Sn  and  PnlLg  ratios  measured  in  the  8-10  Hz  frequency  band.  Each  point  in  these  plots 
represents  the  integrated  gradient  measurement  over  10  km  distance  intervals.  As  shown,  the 
points  nearly  line  up  on  a  linear  trend,  which  we  fit  by  least  squares.  This  same  linear  trend  was  fit 
for  the  other  8  frequency  bands. 


Figures  5(a)  and  (b)  show  the  family  of  curves  which  were  derived  from  this  analysis  for 
Pn/Sn  dindPn/Lg  in  the  nine  filter  bands.  All  frequency  bands  have  the  same  exponentially 
increasing  trend  with  distance,  as  required  by  (12)  and  (13).  The  highest  frequencies,  filter 
passband  above  7  Hz,  seem  to  have  the  smallest  distance  dependence.  The  lowest  frequency  in  the 
Pn/Sn  ratio  (1.5  Hz)  has  the  grestest  frequency  dependence  which  may  be  suspect  because  this 
frequency  band  generally  has  the  lowest  SNR  for  Pn.  In  the  case  of  the  Pn/Lg  ratios,  the  greatest 
frequency  dependence  shows  up  in  the  bands  centered  between  1.5  and  5  Hz. 


Figures  6  (a)  and  (b)  show  the  slopes  of  the  linear  log-ratios  versus  distance  fits  as  a 
function  of  center  frequency  for  Pn/Sn  and  Pn/Lg  ratio,  respectively.  There  is  a  rough  trend  of  the 
slope  decreasing  with  increasing  frequency,  although  there  is  considerable  scatter,  particularly  in 
the  middle  frequency  bands.  Some  of  this  scatter  may  be  due  to  inconsistencies  in  phase  picking 
and  sparsity  of  phase  picks  for  some  distances.  It  should  be  noted  that  only  15  events  were  used 
in  this  preliminary  analysis. 

Although  our  data  fits  produced  somewhat  variable  trends,  they  all  predict  an  exponential, 
or  concave  up,  distance  trend  for  all  frequencies.  These  trends  are  significantly  different  than  the 
trends  from  (5)  and  (6)  based  on  the  Sereno  (1991)  relations,  which  sometimes  predicted  concave 
down  trends,  most  notably  for  the  Pn/Sn  ratios.  We  call  this  method  the  Exponential  Method 
because  it  assumes  that  amplitude  ratios  increase  exponentially  with  distance.  Since  our  curves  are 
based  on  fits  to  actual  data  in  Scandinavia,  an  exponential  distance  trend  best  fits  both  Pn/Sn  and 
Pn/Lg  ratios  in  Scandinavia  over  distance  distance  range  of  200  to  2(XX)  km. 
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Figure  6:  Exponential  attenuation  coefficients  of  the  log  max  rms  Pn/Sn  (a)  and  Pn/Lg  (b)  amplitude  ratios  plotted 
versus  frequency.  Generally,  the  coefficients  are  larger  for  lower  frequency  bands,  although  there  is  considerable 
scatter.  These  coefficients  apply  to  corrections  at  all  distances  in  the  Scandinavian  shield. 


3.3  Distance  Dependence  Correction  Approach 


We  have  implemented  in  ISEIS  a  distance-correction  scheme  for  PIS  amplitude  ratios, 
using  both  approaches  discussed  above.  Figure  7  shows  the  database  indexing  scheme  we  have 
developed  for  accessing  and  applying  path  corrections  for  amplitude-ratio  features.  The  databasing 
scheme  of  ISEIS  is  based  on  the  original  CSS  schema  (Anderson  et  al,  1990)  with  extensions 
described  by  Baumgardt  et  al,  (1991).  The  source  parameters  of  seismic  events  which  have  been 
located  by  IMS  are  stored  in  the  ORIGIN  relation,  indexed  by  a  unique  Orid.  As  shown  in 
Figure  7,  we  assign  events  to  reference  regions,  described  by  the  REFREGION  relation, 
indexed  by  a  unique  Regid.  Reference  events  arc  defined  by  REFEVENT  that  has  the  Orid  of  the 
event  and  the  name  of  the  station,  with  abbreviation  Sta.  Thus,  REFEVENT  assigns  a  specific 
source-receiver  path  to  the  reference  region.  We  add  a  relation,  TECTYPE,  which  gives  a 
tectonic  name  to  the  path  from  the  station  to  the  reference  region.  The  TECTYPE  names  shown 
are  "Shield,"  "Hercynian,"  "Caledonian,"  and  "Cordileran,"  indexed  by  Tectid.  The  tectid  is 
assigned  to  unique  feature  correction  relation,  FEATCORR,  which  defines  a  method  for 
correcting  the  feature  for  distance.  Currently,  we  only  have  one  TECTYPE  instance,  "Shield" 
(for  Scandinavia),  and  two  methods  for  correcting  PnISn  and  PnlLg  amplitude  ratios  for  distance, 
i.e.,  the  Sereno  and  Exponential  method,  described  above.  These  methods  are  algorithms  used 
to  correct  the  amplitude  ratios  to  a  reference  distance. 

This  indexing  scheme  allows  us  to  have  any  number  of  feature-correction  methods 
assigned  to  any  number  of  tectonic  types.  Thus,  in  our  case,  we  have  the  two  methods  discussed 
above  assigned  to  the  "Shield"  tectonic  type,  which  is  Scandinavia.  We  believe  many  feature 
correction  schemes  may  be  generic  for  specific  tectonic  types;  i.e.,  amplimde-ratio  path  correction 
schemes  which  apply  to  the  Scandinavian  shield  will  also  apply  to  other  shield  regions  as  well. 
Future  research  using  data  from  other  regions  of  the  world  will  test  whether  this  is  true.  However, 
our  scheme  allows  new  methods  to  be  added  as  they  are  developed  when  more  data  from  other 
regions  becomes  available. 

Figure  8  shows  an  example  of  the  application  of  the  Exponential  distance  correction 
scheme  to  mesurcments  of  Pn/Sn  amplitude  ratios,  measured  in  the  8  to  10  Hz  band  for  events,  for 
the  seismic  events  near  Novaya  Zemlya  studied  by  Baumgardt  (1993)  and  Fisk  et  al  (1993),  in 
their  analyses  of  the  December  31,  1992  event.  This  event  is  plotted  as  an  asterisk  at  a  distance  of 
about  1 100  km  from  ARCESS  on  the  Pn/Sn  (8-10  Hz)  versus  distance  plot  in  Figure  8.  In  Figure 
8(a),  the  predicted  Exponential  curve  is  plotted  on  the  points,  which  is  the  9  Hz  curve  in  Figure 
5b.  Note  that  on  this  plot,  the  December  31  event  falls  above  the  earthquake  group  and  overlaps 
the  explosion  group  at  a  smaller  distance  from  ARCESS.  Note  that  although  the  curve  seems  to 
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Figure  7:  The  P/S  ratio  correction  indexing  scheme  for  a  station  uses  the  Oracle  database  to  store  tables  of 
reference  events  “indexed”  to  specific  regions  (REFREGIONs)  which  in  turn  point  to  tectonic  types  (TECTYPEs) 
that  classify  the  tectonic  region  containing  the  refregion-  to-  receiver  path.  Tectonic  types  index  relations  that 
point  to  P/S  corrections  (FEATCORRs)  which  are  methods  for  correcting  the  feature  for  distance.  This 
indexing  scheme  allows  us  to  study  a  number  of  different  methods  for  making  corrections. 
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Figure  8:  Examples  of  application  of  the  exponential  distance  corrections  for  Pn/Sn  ratios  measured  at 
ARCESS  (From  Baumgardt,  1993).  (a)  Plot  of  the  distance  correction  curve  for  the  8-10  Hz  band,  (b) 
Same  points  as  (a)  but  with  the  distance  dependence  of  the  curve  removed.  The  correction  moves  the  3 
December  1992  Novaya  Zemlya  event  more  into  the  earthquake  group. 


pass  through  many  of  the  points,  we  are  primarily  interested  in  the  shape  of  this  curve  for 
correction  purposes,  not  its  intercept  with  the  vertical  axis.  Figure  8(b)  shows  the  result  of 
removing  the  exponential  trend  from  the  data  points  so  as  to  shift  the  measurements  to  a  standard 
distance  of  700  km.  When  this  trend  is  removed,  we  find  that  the  December  31  point  now  shifts 
down  and  is  more  consistent  with  the  PnISn  ratios  of  the  earthquakes. 

Figures  9(a)  and  (b)  show  the  same  analysis  using  the  Sereno  curve.  Note  in  Figure  9(a) 
that  the  curve  is  concave  down,  as  we  discussed  above,  although  the  curve  does  seem  to  pass 
through  many  of  the  points.  However,  the  trend  in  the  Sereno  curve  over  this  distance  range  is 
not  signifcantly  different  than  our  Exponential  curve.  Figure  9(b)  shows  the  points  corrected  to 
the  distance  of  700  km  using  the  Sereno  curve.  Again,  the  December  31  point  is  shifted  down  to 
be  more  in  the  earthquake  category,  as  has  also  been  shown  by  Fisk  (1994).  Although  the  shapes 
of  the  Exponential  and  Sereno  correction  curves  are  very  different,  the  resultant  corrections  in 
this  distance  range  have  about  the  same  effects. 

The  scatter  in  the  data  points  in  these  plots  is  too  large  to  make  a  case  for  either  of  the 
Exponential  or  the  Sereno  curves.  This  large  scatter  is  caused  by  the  separation  of  explosions 
and  earthquakes  provided  by  the  high  frequency  Pn/Sn  ratios.  In  fact,  when  the  amplitude  ratios 
are  plotted  on  a  log  scale,  which  we  normally  do,  both  these  curves  are  very  similar  since,  in  log 
space,  they  are  linear  with  about  the  same  slope.  Thus,  when  extrapolating  over  short  distances, 
either  curve  could  be  used  for  distance  correction.  In  the  next  section,  we  show  additional 
examples  of  these  distance  correction  curves  in  log  space  for  features  extracted  from  other  regions. 

We  emphasize  that  both  of  these  sets  of  path  correction  curves  should  be  regarded  as 
preliminary  and  should  only  be  used  with  care.  Sereno's  study  (from  which  the  Sereno 
correction  curves  were  derived,  was  an  assessment  of  regional  detection  capabilities  of  seismic 
networks  for  individual  regional  phases  and  did  not  directly  address  distance  corrections  for 
regional  P/S  amplitude  ratios.  Moreover,  his  study  made  assumptions  about  frequency-dependent 
magnitude  source  scaling.  The  Exponential  approach  does  not  require  source-scaling 
assumptions.  However,  our  Exponential  curves  are  only  based  on  analysis  of  15  events,  and 
more  data  needs  to  be  analyzed  to  refine  the  curves.  Also,  we  relied  entirely  on  the  IMS  phase 
idenfications  in  generating  these  corrections.  We  plan  later  to  verify  the  phase  picks  and  redo  the 
analysis,  if  necessary. 

However,  both  these  sets  of  curves  represent  a  first  step  in  developing  path  corrections  for 
the  amplitude  ratios.  Both  these  studies  clearly  show  that  amplitude  ratios  for  Pn/Sn  and  Pn/Lg 
increase  significantly  with  distances  ranges  over  hundreds  of  km.  Thus,  when  comparing 
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Figure  9:  Same  scatter  plot  as  in  Figure  8  but  showing  the  Pn/Sn  correction  based  on  Sereno  (1991). 

(a)  Plot  of  the  correction  curve  plotted  on  the  scatter  plot.  Note  that  the  curve  trend  is  bending  down 
rather  than  increasing  exponentially,  (b)  After  the  ratios  have  been  corrected  to  a  distance  of  700  km.  As 
in  Figure  8,  after  correction,  the  Pn/Sn  ratio  for  the  31  December  1992  Novaya  Zemlya  event  is  closer 
to  the  earthquakes. 


explosions  and  earthquakes  which  may  be  several  hundreds  of  km  apart,  this  differental  attenuation 
of  Pn,  Sn,  and  Lg  clearly  needs  to  be  corrected  before  the  amplitude  ratios  can  be  used  for 
discrimination. 
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4.0  TRANSPORTABILITY  OF  THE  PnlLg  DISCRIMINANT 


In  the  previous  section,  we  showed  how  to  correct  regional  PIS  amplitude  ratios  for 
distance,  given  that  enough  data  could  be  analyzed  to  generate  a  reliable  set  of  frequency-dependent 
distance  trends  for  all  the  phase  ratios  of  interest.  However,  we  would  expect  that  these  curves 
might  be  quite  different  for  different  tectonic  regions.  For  the  Scandiavian  shield,  attenuation  may 
be  low  for  both  Sn  and  Lg.  However,  for  a  more  tectonically  active  region  of  the  world,  where 
there  is  active  earthquake  and  volcanic  activity,  and  hence,  higher  heat  flows,  we  expect  seismic 
attenuation  in  the  crust  to  be  higher.  Thus,  for  a  given  distance  and  frequency  band,  we  might 
expect  PnlLg  ratios  to  be  larger  in  the  tectonically  active  regions  than  for  the  shield  at  the  same 
distance  since,  for  the  former,  Lg  would  be  more  attenuated.  Moreover,  we  would  expect  the 
distance  dependence  trend  to  be  different  for  shield  and  tectonically  active  regions. 

The  other  problem  is  the  so-called  "Lg  blockages"  which  have  been  observed  in  continental 
areas,  where  tectonic  boundaries  of  different  kinds,  such  as  sedimentary  basins  (Baumgardt, 
1990),  can  either  partially  or  completely  eliminate  the  Lg  energy.  Blockages  are  difficult  to  correct 
by  simple  amplitude-distance  curves  because  they  may  be  very  sudden  and  unpredictable. 
However,  it  may  be  possible  to  find  where  such  blockages  occur  by  exarrxining  crustal  cross 
sections  for  the  source-receiver  propagation  paths  (Baumgardt,  1990).  Once  blockages  are  found 
the  problem  is  how  to  calibrate  Pllg  ratio  discriminants  for  blockages. 

Thus,  the  problem  of  transporting  discriminants  between  different  geographic  regions 
involves  correcting  the  features  for  predictable  effects  due  to  differences  in  the  tectonics  of  the 
different  regions  and  the  more  unpredictable  effects  like  blockages.  It  will  not  be  possible  to 
transport  discriminants  by  just  analyzing  seismic  data  to  calibrate  a  region,  because  in  many 
regions,  seismic  data  may  be  scarce.  In  this  study,  we  have  investigated  using  geological  and 
geophysical  data  now  becoming  available  from  Geographic  Information  Systems  (GIS),  for 
regions  of  high  interest,  such  as  the  Middle  East,  to  correct  for  these  effects. 

4.1  Three  Tectonic  Regions 

To  investigate  the  transportability  of  the  PnlLg  ratio  discriminant  between  different  tectonic 
regions,  we  examine  the  variations  of  discriminants  in  three  different  tectonic  regions,  shown  in 
the  maps  in  Figures  10,  1 1,  and  12.  Each  of  the  maps  show  the  locations  of  events  studied,  and  a 
number  of  the  major  tectonic  features  can  be  seen  on  the  topographic  color  relief  maps.  These 
three  regions  have  very  different  geologic  and  crustal  structures. 
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Figure  10(a)  shows  locations  of  mine  blasts  and  earthquakes  recorded  at  the  ARCESS 
array  in  northern  Scandinavia,  where,  as  shown  in  Figure  10(b),  the  propagation  paths  cross 
primarily  stable,  shield-type  tectonic  regions,  although  the  western  part  of  Norway  would  be 
classified  as  a  Caledonian  Miogeosyncline  (Kvale,  1976).  These  events  have  been  analyzed  by 
ISEIS  in  a  number  of  studies  (Baumgardt,  1993b). 

Figures  ll(a,b)  show  propagation  paths  to  the  GERESS  array  from  mine  blasts, 
earthquakes,  and  rockbursts  in  Germany,  Poland,  and  Switzerland.  The  Vogtland  events  in 
Germany  were  originally  studied  by  Wuster  (1993)  and  in  ISEIS  by  Baumgardt  (1993b).  The  rock 
bursts  and  Swiss  earthquake  swarm  events  appeared  in  the  Ground  Truth  Databases  of  Grant  et  al 
(1993)  and  the  rockbursts  were  studied  in  ISEIS  by  Baumgardt  (1993a).  Figure  11(b)  shows  that 
these  paths'  cross  regions  characterized  as  shield  and  Hercynian  orogens.  The  main  tectonic 
feature  in  the  region  is  the  Tomquist  Zone,  which  separates  the  Hercynian  province  from  the 
western  European  platforms  to  the  north  (Watson,  1976).  This  region  has  been  studied 
extensively  by  refraction  profiles  and  shows  a  pronounced  thickening  of  the  crust  (Guterch  et  al, 
1986). 


Figures  12(a,b)  show  propagation  paths  from  earthquakes  and  nuclear  explosions  in  China 
and  Kazakh  to  the  Chinese  CDSN  station,  Urumchi  (WMQ).  The  tectonic  feature  names  on  the 
topographic  relief  map  are  from  Terman  and  Woo  (1967).  Most  of  the  earthquakes  come  from  the 
data  lists  of  events  previously  studied  by  Bennett  et  al  (1989).  These  include  earthquakes  in  Tien 
Shan,  Mongolia,  and  in  the  northern  part  of  Tibet,  and  Russian  nuclear  explosions  at  the  former 
Soviet  test  site  in  Kazakhstan  (Balapan).  We  have  added  some  additional  more  recent  events,  most 
notably,  the  September  29, 1988  nuclear  explosion  at  Lop  Nor,  which  had  a  magnitude  of  4.7.  To 
our  knowledge,  this  is  the  only  Lop  Nor  nuclear  explosion  for  which  seismic  data  has  been  made 
available  from  WMQ.  These  propagation  paths  cross  tectonically  active  mobile  belts.  There  have 
been  many  studies  of  the  tectonics  and  earthquake  mechanisms  in  this  region  of  China  (Tapponnier 
and  Molnar,  1977;  Molnar  and  Topponnier,  1978;  Ni  and  York,  1978;  Tapponnier  and  Molnar, 
1979)  which  have  shown  that  most  of  the  earthquakes  in  these  regions  have  normal  and  strike-slip 
mechanisms  that  are  consistent  with  the  convergence  of  Eurasia  and  India. 

So,  the  crustal  structures  of  the  three  regions  shown  in  Figures  10,  11,  and  12  progress 
from  simple  shield  stmctures  of  Scandinavia  (Figure  10),  with  little  lateral  heterogeneity  in  crustal 
stmcture,  to  more  active  tectonic  structures  with  much  more  heterogeneous  crustal  structures  of  the 
region  around  GERESS  (Figure  1 1)  and  WMQ  in  China  (Figure  12).  In  this  section,  we  compare 
the  PnlLg  ratio  discriminant  in  these  three  regions  and  investigate  to  what  extent  the  complex  heter- 
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designate  earthquakes  and  circles  designate  blasts,  (b)  Topographic  map  showing  propagation  paths 
from  reference  regions  to  the  station.  Most  of  the  propagation  paths  cross  stable  platform  or  shield 
regions. 


Tectonic  Features  and  Propagation  Paths 
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circles  designate  blasts,  (b)  Topographic  map  showing  propagation  paths  from  reference  regions  to  the 
station.  Propagation  paths  cross  shield  and  tectonic  crustal  structures  and  sedimentary  basins. 


Tectonic  Features  and  Propagation  Paths 
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Figure  12:  These  maps  show  locations  of  earthquakes  and  nuclear  explosions  recorded  at  the  CDSN  station 
WMQ.  (a)  Map  showing  region  names  and  source  types.  Square  symbols  designate  earthquakes  and  circles 
designate  blasts,  (b)  Topographic  map  showing  propagation  paths  from  reference  regions  to  the  station. 
Propagation  paths  cross  a  great  diversity  of  tectonic  boundaries,  fold  belts,  and  sedimentary  basins. 


ogeneous  structures  of  the  tectonically  active  regions  affects  the  values  of  the  PnlLg  ratio  and  must 
be  considered  when  transporting  the  discriminant  from  one  tectonic  region  to  another. 

4.2  Basis  for  the  PnlLe  Amplitude  Discriminant  and  Effects  of  Tectonic 
Structure 

Figure  13  exemplifies  the  basis  for  the  PIS  ratio  discriminant  for  identifying  explosions  and 
earthquakes.  This  figure  compares  an  earthquake  and  explosion  waveform  in  the  Vogtland  region 
of  Germany  recorded  at  one  of  the  elements  of  the  GERESS  array.  The  two  events  are  about  the 
same  distance  (about  190  km)  from  GERESS  and  both  events  are  in  the  same  region  which  means 
that  the  seismic  propagation  paths  to  GERESS  are  about  the  same  for  the  two  events.  On  low 
frequency  or  broadband  traces,  both  earthquakes  and  explosions  produce  sizeable  S  waves, 
although  as  shown  in  Figure  13,  the  amplitude  of  the  S  is  far  larger  than  the  P  for  the  earthquake  as 
compared  with  the  explosion,  where  the  amplitudes  are  nearly  the  same.  This  is  much  more 
pronounced  at  high  frequency  (8-10  Hz)  where  the  explosion  P  amplitude  is  nearly  twice  that  of 
the  S  amplitude  whereas  for  the  earthquake,  the  S  amplitude  is  still  larger  than  the  P  amplitude. 
Also,  the  P  wave  is  impulsive  for  the  blast  and  emergent  for  the  earthquake,  an  observation  also 
made  recently  by  Blandford  (1993)  who  suggested  it  may  be  the  basis  of  a  regional  complexity 
discriminant. 

Figure  14  shows  a  similar  comparison  for  the  Lop  Nor  nuclear  explosion  in  China  and  one 
of  the  nearby  earthquakes  in  Tien  Shan  recorded  at  WMQ.  These  two  events  are  about  the  same 
distance  from  WMQ  that  the  blast  and  earthquake  are  fiom  GERESS  in  Figure  13.  The  Tien  Shan 
earthquake  waveform  in  Figure  14  is  very  similar  to  the  Vogtland  earthquake  waveform  in  Figure 
13.  In  both  cases,  the  S  wave  is  larger  than  the  P  wave  and  the  P  onset  is  emergent.  However, 
we  note  that  the  Lop  Nor  event  has  a  much  larger  P  wave  than  S  wave,  even  at  low  frequency. 
The  nuclear  explosion  in  this  case  seems  to  be  just  a  P  wave  source  and  generates  a  very  weak  S 
wave  in  all  frequency  bands.  It  would  appear,  then,  that  the  nuclear  explosion  in  China  produces 
far  less  shear  wave  energy  than  a  mine  blast  in  Germany.  Either  nuclear  explosions  in  general  do 
not  produce  as  strong  shear  waves  as  mine  blasts  or  something  in  the  propagation  paths  from  Lop 
Nor  to  WMQ  attenuates  the  shear  waves  from  the  Chinese  nuclear  explosion. 

To  compare  nuclear  explosions  in  two  very  different  tectonic  regions.  Figure  15  shows  a 
comparison  of  a  nuclear  explosion  at  Degelen,  recorded  at  WMQ,  and  a  Novaya  Zemlya  nuclear 
explosion,  recorded  at  one  of  the  elements  of  the  ARCESS  array.  The  two  propagation  paths 
cover  about  the  same  distance,  KXX)  km  for  Degelen  to  WMQ  and  1 100  km  for  Novaya  Zemlya  to 
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Figure  13:  Examples  of  sz  waveforms  from  a  Vogtland  mine  blast  and 
earthquake  recorded  at  the  GERESS  array  element  GEC2.  (a)  Mine  blast 
broadband  seismogram  (no  filter),  (b)  earthquake  broadband  seismogram 
(no  filter),  (c)  mine  blast  seismogram,  8-10  Hz  filter,  (d)  earthquake 
seismogram,  8-10  Hz  filter.  At  low  frequency,  both  explosions  and 
earthquakes  generate  strong  S  waves.  At  high  frequency,  explosions  have 
much  less  S  wave  energy  than  earthquakes. 
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Figure  14;  Examples  of  sz  waveforms  from  Lop  Nor  nuclear  explosion 
and  a  Tien  Shan  earthquake  recorded  at  the  CDSN  station  WMQ.  (a)  Nuclear 
explosion  broadband  seismogram  (no  filter),  (b)  earthquake  broadband 
seismogram  (no  filter),  (c)  nuclear  explosion  seismogram,  8-10  Hz  filter,  (d) 
earthquake  seismogram,  8-10  Hz  filter.  At  both  low  and  high  frequency,  the 
earthquake  generates  strong  S  waves  but  the  nuclear  explosion  produces 
weak  S  waves. 
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ARCESS,  but  the  crustal  structures  of  the  paths  are  very  different.  The  propagation  path  cross 
sections  for  the  two  paths,  obtained  from  the  Cornell  crustal  structure  database  for  Eurasia 
(Fielding  et  al,  1992),  are  also  shown  in  Figure  15  and  show  the  variations  in  topography,  depth 
of  sediments,  and  depth  to  the  Moho  along  the  two  paths.  Comparing  the  waveforms  reveals  that 
the  Degelen  nuclear  explosion  has  strong  Lg  waves  at  low  frequency,  but  no  Sn  wave  is  apparent. 
However,  at  high  frequency,  only  a  P  wave  is  apparent.  Evidently,  the  Lg  wave  is  attenuated  at 
high  frequency  and  the  Sn  wave  is  unobserved  at  all  frequencies  in  China.  However,  for  the  path 
from  the  Novaya  Zemlya  nuclear  explosion  to  ARCESS,  which  is  a  shield  path,  there  is  no  Lg 
wave  observed  but  a  strong  Sn  wave  is  observed  at  both  low  and  high  frequency.  The  propagation 
path  cross  section  for  the  Novaya  Zemlya-to- ARCESS  path  shows  that  there  is  a  deep  sedimentary 
basin  in  the  path  which  apparendy  blocks  the  Lg  wave  at  all  frequencies.  However,  Sn  propagates 
efficientiy  in  all  frequency  bands. 

These  comparisons  highlight  the  essence  of  the  discriminant  transportability  problem.  The 
regional  phases  Sn  and  Lg  vary  in  their  observability  in  different  tectonic  regions,  and  also  their 
amplitudes  will  vary  differendy  with  distance  and  frequency,  for  reasons  which  are  not  always  the 
same  in  all  tectonic  regions.  The  crust  in  China  does  not  completely  block  Lg  waves,  although  the 
anelastic  attenuation  of  Lg  in  the  crust  seems  to  be  high  enough  to  eliminate  Lg  at  frequencies 
above  5  Hz.  However,  mantie  attenuation  beneath  China  does  not  appear  to  permit  Sn  propagation 
over  large  distances.  High-frequency  Sn  does  propagate  efficiendy  across  the  Barents  Sea, 
perhaps  because  anelastic  attenuation  is  low  in  the  mantle  beneath  the  region.  However,  Lg  is 
completely  blocked  at  all  fi^uencies,  not  by  anelastic  attenuation,  which  is  low  in  this  region,  but 
by  being  captured  in  the  Barents  Sea  sedimentary  basin,  as  discussed  by  Baumgardt  (1990). 
Transporting  the  Pn/Sn  or  PnlLg  ratio  from  one  tectonic  region  to  another  requires  careful 
accounting  for  these  effects  of  the  propagation-path  structure  on  the  amplitudes  of  the  regional 
shear  phases. 

4.3  Analysis  of  the  Lop  Nor  Nuclear  Explosion  -  A  Test  Case  for  Discriminant 
Transportability 

On  September  29,  1988,  a  nuclear  explosion  test  was  conducted  at  the  Chinese  Lop  Nor 
nuclear  test  site  that  had  a  body-wave  magnitude  of  4.7.  The  closest  station  to  the  Lop  Nor  test  site 
is  the  CDSN  station  Urumchi  (WMQ).  Usually,  this  station  is  turned  off  by  the  Chinese  during 
nuclear  tests  at  Lop  Nor.  However,  for  this  test,  the  station  was  left  on  accidentally.  K  the  nuclear 
explosion  had  not  been  announced  by  the  Chinese,  the  event  might  be  considered  a  special  event 
since  it  was  recorded  at  WMQ  which  should  only  be  on  for  earthquakes. 
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Figure  15:  Comparison  of  propagation  paths  and  waveforms  for  nuclear  explosions  at  Kazakh  recorded  at 
WMQ  and  at  Novaya  Zemlya  recorded  at  ARCESS.  In  China,  Lg  waves  propagate  at  low  frequency  but  Sn 
and  high  frequency  Lg  waves  are  attenuated.  In  Scandinavia,  Sn  waves  propagate  at  all  frequencies  but  Lg 
is  blocked  by  the  Barents  sedimentary  basin.  These  kinds  of  regional  propagation  effects  need  to  be 
understood  and  considered  when  transporting  the  P/S  discriminant  from  one  region  to  another. 


The  characterization  of  this  event  relates  to  the  transportability  of  regional  discriminants 
because  we  have  no  other  nuclear  explosions  at  Lop  Nor  recorded  at  WMQ  with  which  to  compare 
the  event.  The  distance  of  the  event  from  WMQ  is  about  240  km.  The  maps  in  Figure  12  show 
that  most  of  the  events  near  Lop  Nor  are  earthquakes.  The  Russian  nuclear  explosions  at  Kazakh 
are  much  more  distant  (about  1000  km)  and  the  propagation  paths  from  Kazakh  to  WMQ  cross 
many  complex  crustal  blocks.  Thus,  we  can  only  compare  the  event  to  nearby  earthquakes  and 
much  more  distant  nuclear  explosions.  To  characterize  this  event,  in  order  to  identify  it,  assuming 
it  had  not  been  announced  as  a  nuclear  explosion,  there  would  be  two  possible  approaches: 

(1)  Compare  the  PnlLg  ratios  of  the  event,  recorded  at  WMQ,  with  nearby 
earthquakes  recorded  there  and  identify  the  event  as  earthquake  if  it  is 
similar  to  the  earthquakes.  If  not,  the  event  might  be  characterized  as 
anomalous  or  an  explosion.  To  identify  the  event  as  a  nuclear  explosion, 
we  must  compare  it  with  the  more  distant  Russian  nuclear  explosions, 
which  will  require  propagation-path  distance  corrections. 

(2)  Compare  the  PnlLg  ratios  of  the  Lop  Nor  event  with  explosions  and 
earthquakes  recorded  at  comparable  distance  by  some  other  station  in  some 
other  part  of  the  world.  For  example,  there  are  earthquakes  and  mine  blasts 
in  the  ARCESS  and  GERESS  regions,  as  shown  in  Figures  10  and  11, 
which  are  in  the  same  distance  range  as  the  Lop  Nor  nuclear  explosion  is 
from  WMQ.  Comparing  the  Lop  Nor  event  recorded  at  WMQ  with  the 
events  recorded  at  ARCESS  and  GERESS  requires  consideration  of  how 
the  differences  in  the  regional  tectonics  of  the  source  and  propagation  paths 
will  affect  the  discriminants. 

This  scenario  may  often  happen  in  an  actual  CTBT  or  NPT  monitoring  situation,  which 
might  involve  a  "first  test"  in  a  country  which  previously  has  not  admitted  testing  in  a  particular 
region.  The  test  may  be  carried  out  in  a  seismically  active  region  in  hopes  the  test  may  be  confused 
with  an  earthquake.  The  Lop  Nor  test  site  is  located  in  eastern  Tien  Shan  on  the  southern  edge  of 
the  Dzungarian  sedimentary  basin  in  the  vicinity  of  a  number  of  seismically  active  faults 
(Tapponnier  and  Molnar,  1979).  We  first  compare  the  Lop  Nor  event  with  earthquakes  in  this 
region  as  well  as  other  more  distant  earthquakes  in  Tien  Shan,  Tibet,  and  Mongolia,  and  nuclear 
explosions  in  Kazakh.  Then,  we  next  compare  the  Lop  Nor  explosion  with  mine  blasts  and 
earthquakes  recorded  by  GERESS  and  ARCESS  at  comparable  distance  but  in  different  regions. 

Comparison  of  the  Lop  Nor  Explosion  With  Nearby  Events 

In  this  section,  we  discuss  the  ISEIS  analysis  of  the  WMQ  recordings  of  the  Lop  Nor 
nuclear  explosion  and  other  events  around  WMQ,  as  recorded  on  both  the  short-period  (40  Hz 
sampling  rate)  and  broadband  (20  Hz  sampling)  vertical-component  channels  at  WMQ,  referred  to, 
respectively,  as  sz  and  bz.  Wherever  possible,  we  utilize  data  recorded  on  the  40  Hz  sz  channels, 
with  Nyquist  fiequency  of  20  Hz,  which  enables  us  to  analyze  ratios  up  to  the  8-to-16  Hz  band. 
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Figures  16  (a),  (b),  (c),  and  (d)  show  scatter  plots  of  the  PnlLg  ratio,  measured  in  the  8  to 
10  Hz  band  as  a  function  of  region,  magnitude,  ftequency,  and  epicentral  distance,  respectively. 
In  each  plot,  the  Lop  Nor  event  is  plotted  as  a  star  and  is  identified  in  Figure  16(a)  as  the 
CURRENT  EVENT.  Such  plots  are  highly  useful  when  characterizing  special  events  in  order  to 
investigate  all  the  factors  which  could  be  causing  variations  in  the  PnlLg  ratios. 

The  first  plot  in  Figure  16(a)  compares  the  Lop  Nor  event  with  the  other  events  on  the  basis 
of  reference  regions.  The  reference  regions  are  identified  as  follows: 

(1)  Lop_Nor  -  This  includes  the  Lop  Nor  event  itself,  indicated  by  the  plus 
symbol,  and  four  earthquakes,  indicated  by  the  squares.  These  earthquakes 
may  be  associated  with  left-lateral  strike-slip  faulting  in  the  western  Gobi 
desert  (Tapponnier  and  Molnar,  1979). 

(2)  Tien_Shan  -  These  earthquakes  are  located  west  of  Lop  Nor,  in  the  Tien 
Shan  fold  and  fault  belt  on  the  northwestern  part  of  the  Tarim  Stable  belt. 

Most  earthquake  mechanisms  in  this  region  are  consistent  with  thrusting 
associated  with  the  collision  of  the  Indian  plate  with  Eurasia  (Ni,  1978). 

(3)  NEJTibet  -  Three  earthquakes  located  in  the  northeastern  part  of  the 
Tibetan  plateau,  including  one  event  in  the  Tsaidam  sedimentary  basin. 

(4)  SW_Mongolia  -  These  earthquakes  appear  to  be  associated  with  the  NW- 
SE  trending  right-lateral  faults  of  the  Altai  foldbelt  (Tapponnier  and  Molnar, 

1979). 

(5)  NWJTibet  -  These  earthquakes  are  located  in  the  northwestern  part  of 
Tibet,  perhaps  associated  with  the  Altyn  Tagh  left-lateral  strike-slip  fault  (Ni 
and  York,  1978). 

(6)  Balapan  -  These  are  nuclear  explosions  of  the  former  Soviet  Union  in  the 
present-day  Kazakh  Republic.  TTie  "Balapan"  region  includes  events  at  both 
the  Degelen  Mountain  and  Shagan  River  regions  of  the  test  site. 

The  region  plot  in  Figure  16(a)  shows  the  loglO  PnlLg  (8-10  Hz)  ratios  plotted  as  a 
function  of  these  regions  in  order  of  increasing  distance  of  the  regions  from  the  Lop  Nor  event. 
Figure  16(b)  shows  the  same  points  plotted  as  a  function  of  distance  from  WMQ.  Since  the  Lop 
Nor  event  is  also  in  the  Lop  Nor  reference  region,  these  two  plots  are  very  similar.  The 
Scandinavian  Exponential  distance-correction  curve  for  the  8-10  Hz  PnlLg  ratio,  which  was 
derived  in  the  analysis  described  in  Section  2.0,  is  overlayed  on  the  points  in  Figure  16(b).  This 
curve  passes  directly  through  the  earthquakes  in  the  distance  range  of  about  4(X)  to  1000  km,  but  it 
underpredicts  the  nuclear  explosions  and  overpredicts  the  Lop  Nor  earthquakes.  The  two  plots 
show  that  most  of  the  earthquakes  have  log  10  ratios  between  -0.5  to  0.4  over  the  distance  range  of 
400  to  1000  km,  with  no  strong  systematic  distance  trend. 
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may  indicate  systematic  frequency  dependence  of  discrimination.  Large  Pn/Lg  ratios  in  nuclear  blasts 
indicate  blasts  are  primarily  P  sources  at  frequencies  above  6  Hz.  (d)  Plot  of  8-10  Hz  Pn/Lg  ratios  vs 
magnitude  for  Lop  Nor  blast  (star)  and  reference  earthquakes  and  nuclear  explosions.  Magnitude  and 
SNR  plots  can  reveal  discriminant  trends  caused  by  magnitude  variations. 


The  striking  anomalies  in  these  plots  are  the  very  large  PnlLg  values  of  the  Lop  Nor  and 
Kazakh  nuclear  explosions  and  the  small  values  of  the  Lop  Nor  earthquakes.  The  waveforms  for 
these  events  were  shown  in  Figures  14  and  15.  As  pointed  out  earlier,  nuclear  explosions  in 
Kazakh  and  China  appear  to  be  primarily  P  wave  sources  as  recorded  at  WMQ.  Thus,  in  the  8  - 10 
Hz  band,  the  large  values  of  PnlLg  ratio  for  both  Lop  Nor  and  Kazakh  are  due  to  the  fact  that  Lg 
is  very  weak  for  the  Lop  Nor  explosion  and  nonexistent  for  the  Kazakh  explosions.  The  small 
values  for  the  Lop  Nor  earthquakes  seem  to  be  caused  by  strong  Lg  excitation  at  distances  less  than 
about  400  km.  As  shown  in  Figure  16(b),  there  appears  to  be  a  discontinuous  change  in  the  PnlLg 
ratio  at  distances  beyond  about  300  to  400  km.  The  strong  Lg  waves  may  be  due  to  resonances 
and  ground-motion  amplification  produced  by  sedimentary  basins  around  WMQ.  This  point  will 
be  discussed  again  below. 

Figure  16(c)  shows  the  PnlLg  ratio  measurements,  averaged  over  event  type,  as  a  function 
of  the  filter  bandpass  applied  to  the  40  Hz  sampled  sz  data.  Again,  the  star  shows  the 
measurements  for  the  Lop  Nor  nuclear  explosion.  The  error  bars  show  two  times  the  standard 
deviation  of  the  measurements  for  all  the  earthquakes  and  nuclear  explosions.  The  filter  bands 
range  from  0.5-2.5  Hz  on  the  low  end  to  8-16  Hz  on  the  high  end  near  the  20  Hz  Nyquist 
frequency.  This  plot  shows  roughly  increasing  ratios  with  increasing  frequency.  The  greatest 
separation  between  the  Lop  Nor  nuclear  explosion  and  the  earthquakes  is  apparent  in  the  filter 
range  centered  about  the  4  to  6  Hz  band. 

Figure  16(d)  shows  the  scatter  plots  of  the  8-10  Hz  PnlLg  ratios  versus  body- wave 
magnitude,  mb.  This  plot  suggests  a  possible  systematic  trend  of  increasing  ratios  with  magnitude 
for  the  nuclear  explosions  but  no  systematic  trend  for  the  earthquakes.  Also,  Figure  16(c)  shows 
that  the  ratios  for  the  Lop  Nor  nuclear  explosion  decrease  above  the  6  to  8  Hz  band,  and  the 
Kazakh  explosions  increase  above  this  band.  As  we  discussed  above,  nuclear  explosions  in  Asia 
are  primarily  P  sources  at  high-frequencies.  In  fact,  above  8  Hz  there  is  no  Lg  for  the  Kazakh 
explosions  recorded  at  WMQ  (see  Figure  15).  Thus,  the  ratios  in  Figure  16(d)  for  the  nuclear 
explosions  are  actually  PnINoise  ratios,  since  there  appears  to  be  no  F  coda  or  Lg  waves  observed 
from  the  nuclear  explosions  above  8  Hz.  Figure  14  shows  that  the  Lop  Nor  nuclear  explosion 
excites  Lg  waves,  but  they  are  very  weak  compared  with  the  earthquakes.  Therefore,  Figure  15(d) 
shows  that  the  apparent  increase  in  the  ratio  is  caused  by  the  change  in  the  P  amplitude,  relative  to 
noise,  in  proportion  to  the  magnitude  of  the  event 

Finally,  Figure  17  compares  distance-dependent  scatter  plots  of  the  6-8  Hz  loglO  PnlLg 
ratio  recorded  in  by  the  20  Hz  bz  (a)  with  the  40  Hz  channel  (b)  at  WMQ.  Superimposed  on  the 
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Figure  17:  Comparison  of  plots  of  Pn/Lg  ratio  measured  on  the  20  Hz  bz  channel  (a)  and  40  Hz  sz 
channel  (b)  at  WMQ  for  earthquakes  nuclear  explosions  in  Eurasia,  The  astericks  is  the  measurement 
for  the  Lop  Nor  nuclear  explosion.  Note  the  sudden  jump  in  the  ratios  for  the  earthquakes  at  just  beyond 
250  km.  Note  that  the  Scandinavian  distance  correction  curve  almost  fits  the  sz  data  for  the  earthquakes 
after  the  discontinuity. 


plots  are  the  Scandinavian  distance  correction  curves  for  the  6-8  Hz  band  taken  from  the  7  Hz 
curve  in  Figure  5b.  Both  plots  show  very  similar  distance  dependent  trends  at  distances  beyond 
500  km.  At  closer  distances,  the  20  Hz  data  have  somewhat  smaller  PnlLg  ratios.  However,  both 
sets  of  measurements  shows  the  clear  anomalous  change  in  the  Pn/Lg  ratio  beyond  250  km. 

Taken  together.  Figures  16  (a),  (b),  (c),  and  (d)  and  Figure  17  indicate  clearly  that  the  Lop 
Nor  nuclear  explosion  would  be  identified  as  a  pure  compressional  source.  The  earthquakes  in  the 
same  region  as  the  Lop  Nor  test  site  appear  to  be  strong  Lg  sources,  perhaps  due  to  signal 
amplification  in  the  sedimentary  basins  surrounding  WMQ. 

Comparison  With  Events  in  Other  Regions 

Figure  18  shows  a  comparison  of  the  PnJLg  (8-10  Hz)  ratios  for  the  Lop  Nor  event  with 
measurements  for  events  in  the  vicinity  of  GERESS  (a)  and  ARCESS  (b).  The  points  in  Figure  18 
(a)  come  from  a  study  which  showed  clear  discrimination  between  earthquakes  and  explosions  in 
the  Vogtland  region  of  Germany  (Wuster,  1993;  Baumgardt,  1993a)  and  comparison  of  two  such 
events  are  shown  in  Figure  14.  We  also  plot  data  points  for  rockbursts  in  Poland,  first  studied  by 
Baumgardt  (1993a)  which  show  that  these  events  look  like  earthquakes.  The  Lop  Nor  nuclear 
blast  recorded  at  WMQ  has  a  comparable  distance  to  the  Vogtland  mine  blast,  and  the  value  of  the 
ratio  is  similar  to  the  larger  values  of  the  Vogtland  mine  blasts. 

Figure  18(b)  shows  the  same  comparison  with  the  events  in  Scandinavia  recorded  at 
ARCESS,  which  includes  blasts  in  the  Kiruna  mine  of  Sweden,  other  mine  blasts  and  earthquakes 
in  Sweden,  and  earthquakes  in  the  Steigen  region  of  Norway.  Recently,  these  events  were  studied 
by  Baumgardt  (1994)  who  showed  that  the  Kiruna  blasts  have  large  Lg  waves,  comparable  to 
those  generated  by  nearby  earthquakes,  which  suggestes  that  mine  blasts  can  induce  large  shear 
waves,  perhaps  due  to  induced  rock  fracturing  or  spallation  caused  by  ripple  firing.  These  events 
have  ratios  which  are  well  below  those  of  the  Lop  Nor  nuclear  explosion,  recorded  at  WMQ  at 
comparable  distance. 

Thus,  these  two  plots  show  that  the  Pn/Lg  ratios  for  the  Lop  Nor  nuclear  explosion  more 
closely  resemble  the  GERESS  recordings  of  mine  blasts  in  Germany  at  comparable  distance  than 
the  ARCESS  measurements  of  earthquakes  and  blasts  in  Scandinavia.  The  lower  values  of  the 
Pn/Lg  ratio  in  the  northern  pan  of  Scandinavia,  as  measured  at  ARCESS,  may  be  due  to  the  fact 
that  the  propagation  paths  have  less  Lg  attenuation  and  blockage.  As  we  showed  in  Figure  15, 
large  Sn  waves  are  observed  at  high  frequency  from  nuclear  explosions  recorded  in  Scandinavia 
but  not  at  WMQ.  However,  the  Hercynian  Province  of  Germany  appears  to  be  similar  to  the  thrust 
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and  fold  belts  of  China  in  terms  of  the  high  attenuation  of  Lg.  Thus,  we  would  conclude  from  this 
comparison  that  discriminants  are  more  transportable  from  the  GERESS  region  to  China  than  from 
Scandinavia  to  China. 

4.4  Discriminant  Transportability  and  Crustal  Structure 

Crustal  Cross  Sections 

In  this  section,  we  investigate  if  knowledge  of  propagation-path  crustal  stracture  could  be 
utilized  to  transport  discriminants  from  one  tectonic  region  to  another.  This  idea  was  first 
considered  by  Baumgardt  (1990)  who  showed  that  Lg  blockage  could  be  predicted  by  studying 
crustal  cross  sections  of  the  propagation  paths  between  the  source  and  receiver,  like  those  shown 
in  Figure  15.  The  study  revealed  Lg  blockage  in  the  Barents  Sea  sedimentary  basin  between  the 
Novaya  Zemlya  nuclear  explosions  and  the  NORESS  and  ARCESS  arrays.  More  recent  studies 
(Zhang  and  Lay,  1994  a,b;  Zhang  et  al,  1994)  have  shown  correlations  between  crustal  structure 
features  extracted  from  GIS  databases  (Fielding  et  al,  1992)  and  PILg  ratio  measured  for  stations 
over  distances  ranges  of  up  to  20(X)  km.  Strong  correlations  were  observed  between  ratios  in  the 
0.5  to  3.0  Hz  band  and  variations  in  sediment  thickness,  topography,  and  measurements  of  Lg 
attenuation.  However,  because  of  the  large  distances  involved,  only  low  frequencies  could  be 
observed. 

In  this  study,  we  apply  the  method  of  Zhang  et  al  (1994)  to  ISEIS  measurements  of  PnlLg 
ratios  in  frequency  bands  up  to  10  Hz  and  look  for  correlations  with  crustal  structure  features.  In 
order  to  ensure  that  both  Pn  and  Lg  waves  are  above  the  noise  levels  across  this  frequency  band, 
we  used  only  earthquakes  in  this  analysis.  The  crustal  structure  features  were  derived  from  the 
Eurasian  GIS  database  of  Fielding  et  al  (1992).  Our  intent  in  this  study  is  to  determine  if  there  arc 
measurable  features  in  the  propagation-path  crustal  stnicture  which  correlate  with  the  PnlLg  ratio 
features.  If  so,  the  PnlLg  ratio  feature  may  be  "calibrated"  to  adjust  the  ratio  when  transporting  it 
from  one  tectonic  stracture  to  another  very  different  one.  For  the  three  structures  in  Figures  10, 
11,  and  12,  we  seek  crustal  parameters  which  correlate  with  the  PnlLg  ratio  and  can  predict  a 
"baseline"  shift  to  adjust  the  PnlLg  ratio  for  the  Lop  Nor  event  to  allow  it  to  be  compared  with 
explosions  in  recorded  at  GERESS  in  Germany  or  ARCESS  in  Scandinavia. 

Another  difference  between  our  study  and  that  of  Zhang  et  al  (1994)  is  that  we  only  study 
earthquakes  because  we  wish  to  look  for  correlations  between  crustal  stracture  parameters  and 
PnlLg  ratios  at  frequencies  up  to  10  Hz.  Because  the  nuclear  explosions  in  China  appear  to  have 
little  or  no  Lg  energy  in  this  band,  excluding  the  explosions  ensures  that  there  is  sufficient  signal- 
to-noise  ratio  in  the  Lg  at  high  frequencies. 
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Figures  19,  20,  and  21  show  crustal  cross  sections  from  different  "reference  regions," 
defined  by  the  seismic  events  in  the  regions,  to  the  ARCESS,  GERESS,  and  WMQ  stations, 
respectively.  These  cross  sections  indicate  the  variations  in  topography,  depth  to  sediments,  if 
sediments  are  present,  and  depth  to  the  Moho.  The  lines  on  the  topographic  maps  between  the 
stations  and  the  specified  regions  are  shown  in  Figures  10,  11,  and  12.  For  this  analysis,  we 
categorize  our  events  into  geographic  clusters  using  the  ISEIS  reference  region  concept, 
illustrated  in  Figure  7  and  discussed  in  Section  3.3.  The  reference  regions  have  centroid 
locations  determined  by  the  average  location  of  clusters  of  events  in  the  cluster  and  stored  in  the 
Oracle  database  in  the  REFREGION  relation  indexed  by  a  unique  regid  number  and  a  textual 
regname.  Specific  seismic  events,  stored  in  the  database  in  the  ORIGIN  relation  indexed  by 
orids,  are  assigned  to  the  reference  region  through  the  REFEVENT  relation,  which  associates 
the  arid  and  the  regid  and  sta.  (see  Section  3.3  and  Figure  7).  In  our  analysis  below,  we  examine 
the  effects  of  propagation  paths  to  reference  region  centroids  on  the  average  PnlLg  ratio  for  the 
events  in  the  region. 

In  Figure  19,  the  three  propagation  paths  are  from  ARCESS  to  the  average  location  of  the 
Kiruna  mine  blasts  reference  region  (Kiruna),  earthquakes  and  explosions  in  northern  Sweden 
reference  region  (NJSweden),  and  the  average  location  of  the  Steigen  earthquake  reference 
region  (Steigen).  Generally,  these  cross  sections  indicate  that  the  crustal  cross  sections  to  all 
these  reference  regions  have  uniform  thickness  of  about  40  km,  with  little  topography  and  no 
sediments. 

The  reference  regions  in  Figures  1 1  and  20  for  GERESS  are  Switzerland,  defined  by  an 
earthquake  swarm  in  Switzerland,  Upper _Silesia,  defined  by  rockbursts  in  southeastern  Poland, 
Lubin,  defined  by  rockbursts  in  south-central  Poland,  and  Vogtland,  defined  by  the  Vogtland 
earthquakes  and  blasts.  These  cross  sections  exhibit  significantly  more  variation  in  crustal 
structure,  as  compared  to  the  Scandinavian  region  around  ARCESS,  particularly  in  the  crustal 
thickness.  The  crust  is  generally  about  30  to  35  km  for  most  of  these  paths  to  GERESS,  except 
for  the  paths  which  cross  the  Swiss  Alps,  where  the  crust  thickens  to  50  km. 

Finally,  Figures  12  and  21  show  the  propagation  paths  and  cross  sections  between  WMQ 
and  the  reference  regions  in  China,  discussed  above  (Lop_Nor,  NW_Tibet,  NEJTibet, 
WJTienshan,  and  Mongolia).  The  former  Soviet  testsite,  called  Balapan,  is  also  shown, 
although  in  this  study,  we  exclude  nuclear  explosions).  These  regions  show  the  greatest 
complexity  in  surface  topography  and  sediment  thickness,  because  the  propagation  paths  in  Figure 
22  cross  many  foldbelts  and  sedimentary  basins.  Most  notable  are  the  sedimentary  basins  located 
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cross  a  small  sedimentary  basin,  where  sediments  thickness  reach  3  to  4  km. 
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result  in  anomalously  low  Pn/Lg  ratios. 


around  WMQ,  which  are  the  fault-bounded  Dzungarian  basin  to  the  north  and  Tarim  basin  to  the 
south. 


It  should  be  noted  that  the  accuracy  of  these  cross  sections  is  variable.  Although  the 
Cornell  group  has  tried  to  assure  the  reliability  of  the  databases  placed  on  their  Internet  server,  the 
accuracy  of  the  database  is  not  always  certain.  For  example,  the  digitized  sediment  thickness 
maps,  from  which  the  sediment  thickness  cross  sections  were  derived,  are  gridded  to  5'  by  5' 
squares.  For  a  particular  cross  section,  this  yields  approximately  9  km  accuracy  in  the  distance 
position  of  the  sediment  thickness.  We  note  in  particular  that  the  indication  of  thick  sediments 
under  WMQ  is  uncertain,  since  WMQ  appears  to  lie  on  resistant  sedimentary  rocks  in  the  foothills 
north  of  the  Dzungarian  basin  (Francis  Wu,  personal  communication).  However,  the  Eurasian 
database  indicates  that  it  lies  on  up  to  8  km  of  sediments  above  the  crystalline  basement,  because  of 
the  low  resolution  of  the  sediment  thickness  database.  Because  the  position  of  these  sediments 
may  be  off  by  as  much  as  9  km,  it  would  not  be  accurate  to  place  WMQ  inside  a  sedimentary  basin 
itself,  although  it  definitely  lies  within  8  or  9  km  of  the  edge  of  the  Dzungarian  basin. 

With  the  caveat  that  the  present  accuracy  of  the  geological  databases  is  limited,  we  present 
our  approach  to  using  crustal  parameter  correlations  to  calibrating  the  PnlLg  ratio  discriminant  for 
different  tectonic  regions.  The  accuracy  of  the  actual  calibration  in  a  region  should  improve  as  the 
accuracy  of  the  GIS  databases  improve.  At  least  this  approach  may  provide  a  first  order  estimate 
of  the  calibration  for  regional  tectonics. 

Crustal  Structure  Parameters 

Following  Zhang  et  al  (1994),  we  define  the  following  crustal  structure  parameters  derived 
from  the  propagation-path  crustal  cross  sections: 

(1)  £(A;)  -  The  elevation,  in  meters,  at  the  distance  A/  from  the  station. 

(2)  CT(A. )  -  The  crustal  thickness,  in  kilometers,  at  the  distance  A/  from  the  station 

(3)  57’(  A; )  -  The  sediment  thickness,  in  kilometers,  at  the  distance  A/,  from  the  station 

These  parameters  are  then  averaged  over  each  distance  point  along  the  cross  section  for  the  entire 
path. 


Baumgardt  (1990)  suggested  that  it  was  the  change  in  crustal  parameter,  in  particular, 
changes  in  the  sediment  thickness  above  the  basement,  more  than  the  mean  value  of  the  parameter 
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itself,  which  might  cause  blockage  of  Lg  waves  in  continental  areas.  Zhang  et  al  (1994) 
considered  parameterizations  of  the  maximum  and  minimum  values  of  elevation  and  crustal 
thickness  and  maximum  sediment  thickness,  in  addition  to  the  mean  values.  However,  the 
maximum  and  minimum  values  may  be  redundant  with  the  mean  values.  In  this  study,  we 
investigate  a  different  parameterization  which  is  a  more  direct  measure  in  the  change  in  these 
features  with  distance,  called  the  average  upper-quartile  of  the  gradient.  To  compute  this 
parameter  for  each  path,  we  first  compute  the  gradient  of  the  parameter  at  each  distance  point  along 
the  path  as  central  difference  derivatives,  as  follows: 


(4) 


V[£(A,)] 


£(A,,,)-£(A,.,) 

^<■+1  ~  ^i-l 


-  The  elevation  gradient  at  distance  point.  A/,  in 


units  of  meters/kilometers. 


(5) 


V[CT(A  )]  =  -  The  crustal  thickness  gradient  at  distance.  A/ 

in  units  of  kilometers/kilometer. 


(6) 


V[5r(A  )]  =  i)  .  sediment  thickness  gradient  at  distance,  Aj 


in  units  of  kilometers/kilometer. 


For  each  distance  point  along  the  path,  these  central  difference  gradient  values  are  computed. 
Rather  than  average  these,  we  determine  the  upper  quartile  of  these  derivatives  and  average  these 
values.  Thus,  the  estimated  gradient  is  an  estimate  of  the  largest  gradients  observed  along  the  path. 
Averaging  the  upper  quartile,  rather  than  taking  the  maximum  gradient,  helps  to  avoid  spuriously 
large  values  due  to  the  limited  resolution  along  the  path. 

Correlation  Analysis 

In  the  approach  of  Zhang  et  al  (1994),  crustal  cross  sections  were  determined  for  all  paths 
between  each  source  and  receiver,  and  the  parameters  for  each  path  were  computed.  As  discussed 
above,  we  compute  the  crustal  parameters  for  the  propagation  path  to  each  reference  region  and 
examine  the  correlation  between  the  mean  values  of  the  Pn/Lg  ratio  for  the  region  and  the  crustal 
parameters  for  the  path  between  the  station  and  the  reference  region.  Each  of  the  parameters 
defined  above  are  computed  for  each  of  the  crustal  cross  sections  shown  in  Figures  19,  20,  and  21 
and  ISEIS  used  to  measure  and  average  the  PnlLg  ratios  for  the  reference  regions. 

Figures  22, 23,  and  24  show  correlation  plots  of  the  average  PnlLg  ratio,  measured  in  four 
frequency  bands,  0.5-2.5  Hz,  2-4  Hz,  5-7  Hz,  and  8-10  Hz,  as  a  function  of  the  mean  values  of 
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Figure  22:  Correlation  scatter  plots  of  the  log  10  Pn/Lg  amplitude  ratios 
versus  mean  elevation  along  the  propagation  path. 
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Figure  23:  Correlation  scatter  plots  of  the  log  10  Pn/Lg  amplitude  ratios 
versus  crustal  thickness  along  the  propagation  path. 
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Figure  24:  Correlation  scatter  plots  of  the  log  10  Pn/Lg  amplitude  ratios 
versus  mean  sediment  thickness  along  the  propagation  path. 
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elevation,  crustal  thickness,  and  sediment  thickness,  respectively.  Plots  of  the  average  of  the 
upper  quartile  of  the  gradients  of  these  parameters  in  the  same  four  frequency  bands  are  plotted  in 
Figures  25, 26,  and  27.  The  error  bars  on  each  point  are  the  standard  deviations  about  the  means 
of  the  PnlLg  ratio  measurements  for  each  reference  region.  For  each  plot,  a  straight  line, 
determined  by  linear  regression,  is  plotted  through  the  points  and  the  slope  and  standard  errors  of 
the  regressions  are  shown  on  each  plot. 

In  all  these  analyses,  we  find  the  best  correlations  with  the  high-frequency  ratios  (5-7  Hz 
and  8-10  Hz)  both  in  terms  of  size  of  the  slope  and  minimum  error.  The  smallest  standard- 
deviation  error  bars  are  found  for  the  highest  frequencies,  with  the  exception  of  the  Lop  Nor 
region,  where  there  are  large  variations  in  the  Pn/Lg  ratio,  as  we  have  discussed  above.  We 
discuss  this  error  bar  in  more  detail  below.  Most  of  the  error  bars  indicate  a  scatter  of  on  the  order 
of  -t--  0.5  log  PnlLg  ratio  units,  which  is  comparable  to  the  scatter  found  in  the  Zhang  et  al  (1994) 
study. 

At  high  frequency,  the  best  correlated  parameters  are  mean  elevation,  mean  sediment 
thickness,  and  their  upper  quartile  gradients  of  the  parameters.  The  crustal  thickness  parameter 
shows  very  little,  if  any,  correlation  with  PnlLg  ratios.  Baumgardt  (1990)  argued  that  Lg 
attenuation  and  blockage  is  primarily  dominated  by  upper  crustal  effects,  like  the  upper-crustal 
velocity  variations  caused  by  sedimentary  basins  and  elevation  changes.  Depth  to  the  Moho 
appears  to  have  very  little  effect  on  the  variations  in  PnlLg  ratios. 

To  compare  with  the  Zhang  et  al  (1994)  study,  the  0.5  to  2.5  Hz  PnlLg  ratios  are  closest  to 
the  0.5  to  3.0  Hz  PILg  ratios  of  their  study.  Figure  22  agrees  with  the  Zhang  et  al  (1994)  study  in 
that  the  mean  elevation  has  a  slightly  negative  conelation  with  the  log  PnlLg  ratio.  However,  at  the 
higher  frequencies  we  have  examined  in  this  study,  the  correlations  become  positively  sloped. 
Also,  our  study  agrees  with  the  Zhang  et  al  (1994)  study  in  that  the  best  correlations  are  found  with 
sediment  thickness  and  elevation  parameter  and  the  worst  correlation  is  with  the  crustal  thickness. 

We  conclude  that  the  best  correlation  appears  to  be  with  sediment  thickness  variations. 
Figure  28  shows  plots  of  the  loglO  PnlLg  ratio  in  the  8-to-16  Hz  frequency  band  plotted  versus  (a) 
the  mean  sediment  thickness  and  (b)  the  gradient  in  the  sediment  thickness.  Strong  correlations  are 
apparent,  as  indicated  by  the  large  positive  slopes  and  the  small  errors  relative  to  the  size  of  the 
slopes.  Note  that  these  regressions  do  not  include  the  Lop  Nor  point,  which  we  regard  as 
anomalous,  as  shown  in  Figure  29.  This  large  variance  may  be  due  to  large  variations  in  Lg 
excitation  due  to  resonances  in  the  sedimentary  basins  near  WMQ.  We  would  expect  this  to  occur 
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Ave.  Gradient  Elevation  (m/km) 


Figure  25:  Correlation  scatter  plots  of  the  loglO  Pn/Lg  amplitude  ratios 
versus  average  upper  quartile  of  the  elevation  gradient. 
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Figure  26:  Correlation  scatter  plots  of  the  log  10  Pn/Lg  amplitude  ratios 
versus  average  upper  quartile  of  the  crustal  thickness  gradient  along  the 
propagation  path. 
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Figure  27:  Correlation  scatter  plots  of  the  log  10  Pn/Lg  amplitude  ratios 
versus  average  upper  quartile  of  the  sediment  thickness  gradient  along  the 
propagation  path. 
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Figure  28:  These  plots  show  that  correlation  of  8-16  Hz  PnA-g  ratios  with  average  crustal  thickness 
(a)  and  the  average  of  the  upper  quartiie  of  the  sediment  thickness  gradients  (b).  Correlation 
regression  does  not  include  the  Lop  Nor  point.  Anomalous  values  for  the  Lop  Nor  earthquakes  may 
result  from  strong  Lg  amplitudes  due  to  sedimentary  basin  resonances. 


primarily  in  situations  where  either,  or  both  the  source  and  receiver,  are  in  the  same  sedimentary 
basin.  If  the  Lg  must  cross  an  entire  basin,  there  may  be  partial  blockage,  of  the  Lg  energy. 

These  plots,  show  that  the  variations  in  PnlLg  ratios  consistently  depend  on  the  tectonic 
type  for  the  source  region  and  propagation  path.  For  example,  looking  at  Figure  28  for  the  highest 
frequencies,  the  smallest  ratios  (i.e.,  largest  Lgs)  are  found  for  propagation  paths  in  the  shield 
regions  (Northern  Sweden,  Steigen,  Vogtland)  and  the  highest  ratios  (i.e.,  weakest  Lgs) 
correspond  to  the  orogenically  active  areas  (Tien  Shan,  Tibet,  Mongolia).  The  Hercynian  orogenic 
regions  of  Europe  (Switzerland,  Upper  Silesia)  appear  to  be  very  similar  with  certain  regions  of 
China  (NE  Tibet,  Mongolia).  There  is  roughly  a  shift  of  about  0.8  log  units  between  the  shield 
and  the  orogenic  region.  This  shift  seems  to  correlate  with  increasing  sediment  thickness  and 
thickness  variations. 

This  approximate  correlation  suggests,  but  does  not  prove,  a  causal  relation  between 
sedimentary  basin  thicknesses  and  variations  in  thickness  with  the  relative  PnJLg  exciation. 
However,  this  study,  along  with  that  of  Zhang  et  al  (1994),  strongly  suggests  a  correlation  that 
may  be  predictive  of  the  variations  in  the  Pn/Lg  ratio  when  one  moves  from  a  shield  region  to  a 
complex  tectonic  region,  irrespective  of  how  it  is  produced.  Thus,  we  would  shift  the  values  of 
discrimination  thresholds  by  about  0.5  to  0.8  when  transporting  discriminants  from  the  shield 
region  to  the  orogenic  region.  This  would  not  work  for  situations  like  Lop  Nor,  which  we 
consider  to  be  anomalous.  However,  examination  of  crustal  cross  sections  would  indicate  if  such 
conditions  like  those  at  Lop  Nor  exist.  This  baseline  shift  would  apply  for  all  distances  contained 
within  the  region.  Once  inside  the  region,  corrections  for  distance  would  be  made  using 
methodologies  like  those  discussed  in  the  last  section. 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Summary 

This  study  has  shown  that  the  differential  effects  on  P  and  S  waves  of  structure  along  the 
propagation  path  from  source  to  receiver  on  P  and  S  waves  dominate  the  transportability  of  the 
regional  PIS  ratio  discriminant.  First,  the  amplitudes  of  Pn,  Sn,  and  Lg  attenuate  at  different  rates; 
shear-wave  attenuation  is  greater  than  compressional-wave  attenuation,  even  for  paths  in  stable 
shield  regions.  Second,  crustal-structure  anomalies,  such  as  sedimentary  basins,  prirmarily  block 
Lg  waves,  which  propagate  mainly  in  the  upper  crust,  whereas  Pn  and  Sn  propagate  through  the 
mantle  beneath  these  anomalies  and  are  unaffected.  Pn  and  Sn  waves  propagate  mainly  in  the  upper 
mantle  lid  above  the  mantle  low  velocity  zone.  However,  Sn  appears  to  be  strongly  attenuated  in 
tectonically  active  regions,  perhaps  due  to  the  low  Q  low-velocity  20  zone  which  is  more 
pronounced  in  these  regions. 

Transportability  of  the  regional  P/S  ratio  discriminant  will  be  of  great  importance  in  an 
operational  CTBT  monitoring  system,  because  for  most  parts  of  the  world  that  must  be  monitored, 
there  has  been  no  nuclear  weapons  testing.  Identifying  the  "first  test"  requires  using  training  data 
from  parts  of  the  world  where  nuclear  testing  has  occurred.  Thus,  differences  in  the  geology  of 
the  different  regions  wherein  seismic  events  are  being  compared,  in  terms  of  the  effects  of  these 
differences  on  on  P  and  S  propagation,  need  to  be  considered  when  using  the  PIS  ratio  for  event 
identification. 

5.2  Conclusions 

In  this  study,  we  have  addressed  the  transportability  issue  by  comparing  amplitude  ratio 
measurements  for  earthquakes  and  explosions  using  data  from  three  verydifferent  tectonic  regions. 
Based  on  these  analyses,  we  arrive  at  the  following  overall  conclusions: 

•  Distance  correction  curves  can  be  derived  for  continental  areas,  assuming 
enough  events  recorded  at  two  or  more  stations  can  found  that  cover  a  wide 
distance  range  and  can  be  used  to  correct  Pn/Sn  and  Pn/Lg  ratios  to  a 
standard  distance.  We  find  exponential  models  seem  to  best  fit  the  observed 
data  for  mine  blasts  and  earthquakes  in  Scandinavia.  We  have  assumed 
geometric  spreading  to  be  frequency  independent  and  the  same  for  P  and  S 
waves. 

•  Our  analysis  of  the  amplitude  ratio  corrections  from  the  Sereno  (1991) 
detection  study,  used  by  Fisk  (1994)  to  correct  PnISn  ratios  of  the  Novaya 
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Zemlya  event  for  distance,  has  a  "bending  down"  type  distance  trend  for 
PniSn  ratios  at  frequencies  below  14  Hz.  However,  the  PrULg  ratio  curves 
have  a  cubic  dependence  with  distance.  Both  of  these  trends  to  not  seem  to 
be  consistent  with  our  observations.  However,  with  the  limited  data 
currently  available  and  the  large  scatter  in  the  obsovations,  it  is  not  possible 
to  determine  if  either  of  the  correction  curves  is  best.  Thus,  we  have 
incorporated  both  distance-correction  schemes  in  ISEIS. 

•  In  tectonically  active  regions,  such  as  China,  nuclear  explosions  recorded  at 
high  frequencies  at  all  distances  are  fundementally  P-wave  sources.  Lg 
energy  is  strongly  attenuated  at  high  frequency  and  Sn  seems  not  to 
propagate  at  all.  In  shields,  however,  Lg  can  be  blocked  in  the  crust  by 
tectonic  features  like  sedimentary  basins  but  high  frequency  Sn  phases  can 
propagate  to  large  distances. 

•  The  Lop  Nor  nuclear  explosion,  recorded  at  WMQ,  can  be  clearly 
discriminated  from  nearby  earthquakes  on  the  basis  of  its  high  PnlLg  ratio. 
Earthquakes  near  WMQ  excite  very  strong  Lg  waves  and  thus  produce  very 
low  PnlLg  ratios  which  appear  anomalous  compared  to  more  distant 
earthquakes.  These  events  all  occurred  well  inside  a  sedimentary  basin  and 
WMQ  is  near  the  edge  of  a  sedimentary  basin.  Thus,  the  largeLg  waves 
may  have  been  generated  by  resonances  in  the  soft  sediments  of  the  basin. 

The  fact  that  the  Lop  Nor  nuclear  explosion  is  also  in  a  sedimentary  basin, 
but  produced  weak  Lg  waves,  shows  further  how  well  the  PnlLg  ratio 
discriminant  works. 

•  Examination  of  crustal  structures  for  the  three  regions  reveals  that  PnlLg 
ratio  correlates  with  crustal  structure  parameters.  Like  the  study  of  Zhang  et 
al  (1994),  we  find  the  greatest  correlation  with  elevation  and  sediment 
thickness  and  their  gradients.  Moho  depth  variations  have  little  effect  on 
variations  in  PnlLg  ratios. 

•  Comparing  the  PnlLg  ratios  of  the  Lop  Nor  nuclear  explosion  with  those  of 
events  near  ARCESS  and  GERESS,  we  find  that  the  Lop  Nor  values  are 
similar  to  those  of  mine  blasts  in  Vogtland  recorded  at  GERESS.  ARCESS 
recordings  of  mine  blasts  in  Scandinavia  have  much  smaller  PnlLg  ratios 
than  the  Lop  Nor  nuclear  explosion.  Since  ARCESS  is  in  a  shield  region, 
we  attribute  the  lower  PnlLg  ratios  to  the  fact  that  Lg  is  less  attenuated 
relative  to  Pn  in  the  shield. 

5.3  Recommendations 

Further  data  collection  and  research  will  be  required  to  fully  assess  the  methodologies  for 
transporting  regional  discriminants.  We  have  only  studied  three  different  tectonic  regions. 
Analysis  of  data  from  other  tectonic  regions,  including  regions  analogous  to  the  Scandinavian 
shield  and  tectonic  orogens  we  have  studied,  need  to  be  analyzed.  Therefore,  we  recommend  that 
further  research  with  discriminant  features  be  conducted  along  the  following  lines: 
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•  Further  refine  the  distance  correction  curves  for  PnISn  and  PnlLg  ratios  for 
Scandinavia  by  analyzing  more  events.  As  for  China,  the  CDSN  station 
coverage  was  too  sparse  for  us  to  apply  the  same  analysis  there  because  we 
were  unable  to  find  enough  earthquakes  in  China  recorded  at  two  or  more 
CDSN  stations  to  derive  Ae  same  curves.  It  may  be  possible  to  accomplish 
this  with  more  stations  in  the  region. 

•  The  method  of  correlating  PIS  ratios  with  crustal  structure  parameters  is 
promising.  However,  additional  crustal  structures  would  be  useful. 
Geological  and  crustal  structure  information  is  now  becoming  available  for 
the  Middle  East  When  seismic  data  becomes  available  for  this  region  in  the 
future,  similar  studies  there  might  be  possible. 


As  for  tranporting  discriminants,  we  believe  the  method  of  comparing 
events  in  different  regions  can  be  a  very  effective  method  of  characterizing 
events  in  new  regions,  where  there  are  limited  training  events.  Perhaps 
very  simple  "baseline  corrections,"  which  do  not  seem  to  be  very  large, 
would  be  sufficient  for  correcting  ratios  for  different  regions.  In  general, 
comparing  PnlLg  ratios  of  shield  regions  of  Scandinavia  with  those  of 
tectonic  active  regions  of  China  would  require  the  latter  ratios  to  be  shifted 
down  by  as  much  as  0.5  to  0.8  log  units  relative  to  the  shield  values.  These 
baseline  shifts  may  be  derived  for  different  tectonic  types  based  on 
extrapolations  of  the  correlations  of  PnlLg  ratio  with  crustal  parameters. 
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